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Groups of Order 2™ Which Contain a Relatively Large 
Number of Operators of Order Two. 


By G. A. MILLER. 


§1. Introduction 


When exactly one-half of the operators of a group are of order 2 the 
group is obtained by extending an abelian group of odd order by means of 
an operator of order 2 which transforms each operator of this abelian group 
into its inverse, and if a group is constructed in this manner then exactly 
one-half of its operators are of order 2. If more than one-half of the opera- 
tors of a group are of order 2 this excess must always be a number of the 
form 2¢—1 and there is an infinite system of such groups for every positive 
integral value of a.* 

From the fact that the number which expresses the excess of the operators 
of order 2 over half the order of the group is always of the form 2*—1 at 
results directly that in a group of order 2” the number of operators whose 
orders exceed 2 is always equal to the order of the group multiplied by a 


and there are groups in which this product gives 


Qn 
number of the form — 


the number of the operators whose orders exceed 2 for every positive integral 
value of n. The smallest relative number of operators whose order exceeds 2 
in a group of order 2” is therefore one-fourth of the order of the group in 
case there is at least one such operator. 

Let G be any group in which more than one-half of the operators are 
of order 2. When G is abelian it is evidently of type (1,1,1,--- ) and 
hence we shall assume in what follows that G is non-abelian. As more than 
one-half of the operators of G are of order 2 these operators generate @ and 
hence some of the operators of order 2 contained @ must be non-invariant. 
If we let H, represent the subgroup composed of all the operators of G which 
are commutative with one such non-invariant operator s, then more than one- 
half of the operators of H, must be of order 2. 

If H, is non-abelian it contains a non-invariant operator s, of order 2, 
and we may represent by H, the subgroup of H, composed of all its operators 
which are commutative with s,. The subgroup H, includes the central of G 


* Bulletin of the American Mathematical Society, Vol. XXV (1919), p. 33). 
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and H, includes this central as well as s,. If H, is non-abelian we can con- 
tinue this process until we arrive at an abelian group Hm which includes the 
central of G and s,, 8,,°°** , Sm. The group Hm is evidently of type 
1,1,1,-.--: ). In the article to which we referred it was proved that m = 1 
whenever the order of G is divisible by some odd prime number but it is 
possible to construct groups of order 2” in which m exceeds any given 
integer. 

All the operators of G can be uniquely represented in the following form 


where it may be assumed without loss of generality that each of the operators 
t., ts,° °° , tr is of order 2. It is known that whenever the order of @ is 
divisible by an odd prime number then exactly one-half of the operators in 
each of the co-sets Hmta, 2 =6=r, are of order 2 and that in no case can 
more than one of these co-sets involve more operators whose orders exceed 2 
than operators of order 2. In the following section we shall prove that when 


there is one such co-set and 7 > 2 then exactly one-fourth of its operators are 


of order two. 


§ 2. Less than one-half of the operators in one co-set are of order two 


The co-set in which less than one-half of the operators are of order 2 is 
characterized by the fact that it is composed of all the operators of G which 
are commutative with less than one-half of the operators of Hm. In parti- 
cular, Hm contains an operator s, which transforms an operator ¢, of this 
co-set into its inverse and hence if H, is defined as above it does not include 
any operator of this co-set. Every operator of order 2 found in H, but not 
in H,», r > 2, is either commutative with ¢, or transforms ¢, into its inverse, 
for if the product of ¢, and this operator of order 2 has an order which exceeds 
2 then this product is transformed into its inverse by s, and hence the said 
operator of order 2 and ¢, are commutative. 

The operators of order 2 contained in H, but not in H» must either 
generate H, or an abelian subgroup of type (1,1,1, - ~:~: ) and index 2 under 
H,. The former of these alternatives is impossible since ¢, cannot be com- 
mutative with one-half of the operators of Hm. Otherwise one-half of the 
operators of the co-set containing ¢, would be of order 2. It therefore results 
that when r > 2 exactly one-fourth of the operators of H» are commutative 
with ¢, and that the central of H, is a subgroup of index 2 under Hm. As 
exactly one-fourth of the operators of the co-set containing ¢, are of order 2 
we have proved the following theorem. Whenever one of the co-sets 
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Hyt,.,' ++, Hmtr, r > 2, contains more operators of order exceeding 2 than 
of order 2 then exactly one-fourth of the operators of this co-set are of order 2. 

It may be assumed without loss of generality that Hmt, is the co-set in 
which less than one-half of the operators are of order 2, and that H, is com- 
posed of the co-sets Hmt,, ++: , Hmts in addition to Hm. When 82 the 
subgroup H» is invariant under H, because it is of index 2. If 8 > 2 this 
subgroup is still invariant since it is generated by all the operators of order 2 
found in H, and having more than two conjugates under H,. It has there- 
fore been proved that H, is invariant under H, independently of the value 
of 6. 

It is now easy to prove that H» is also invariant under G. In fact, the 
operator s, could have been so selected that it would be commutative with all 
the operators of any one of the given co-sets except the co-set Hmt,, and hence 
H,, is invariant under each of these co-sets. It therefore results that when- 
ever more than one-half of the operators in one of the co-sets Hmta,2 SaSr 
have orders which exceed 2 then Hm 1s invariant under G. 

The quotient group G/Hm involves only operators of order 2 besides 
the identity and hence this quotient group must be abelian and of type 
(1,1,1,--- ). The operators of Hm which are commutative with ¢, con- 
stitute the central of G and include the commutators of G which arise from 
operators in H». The former of the theorems involved in the last sentence 
results directly from the fact that these operators are also commutative with 
every operator in the co-sets Hmta, 2 a> 8, while the latter is a consequence 
of the fact that each of the operators ¢,,--- , ¢- transforms H» into an auto- 
morphism of order 2. If the commutator of order 2 arising from 
ta, 2 senr—t, would not be commutative with ¢, then tat, could not 
transform H» according to an automorphism of order 2. 

Since all the commutators arising from operators of Hm are found in 
the central of G and since each of the operators ¢,, - - - ty_, is commutative 
with exactly one-half the operators of H» it results that all the commutators 
with respect to Hm arising from these operators must also arise from ¢,. As 
only three such commutators of order 2 arise from ft, it results that § = 2 
and r= 4. The central of G may have any order of the form 2”, m’” >1 
and there is evidently one and only one such group for every possible value 
of m’’. Hence the following theorem has been established. There is one 
and only one group of order 2%, a > 5, such that one of the co-sets with respect 
to H» involves more operators whose orders exceed 2 than operators of order 2. 
This group contains an abelian subgroup of index 4 and of type (1,1,1,-- -°) 
and a central of index 16. The central quotient group is abelian and of type 
(1,1,1,--+ + ), and the commutator subgroup is of order 4. 


| 
| 
| 


MittER: Groups of Order 2” which Contain 


When a> 6 this group is evidently the direct product of the group of 
order 64 which belongs to this system and the abelian group of order 2*° and 
of type (1,1,1,--.- ). In the form in which G presented itself above m = 2 
but it is possible to select for s, an operator of Hm so as to make m=—1. If 
this is done the order of Hm remains the same since one of the corresponding 
co-sets must have less than one-half of its of order 2. When half of the 
operators of each co-set must be of order 2 it is clear that the order of Hm 
is independent of the choice of the set of operators s,, 8.,° °° Sm. Hence it 
results from the property of the groups just considered that the order of Hm 
is an invariant of every group in which at least one-half of the operators are 
of order 2 and r > 2. 

It should also be noted that when 7 > 2 and G has a co-set, with respect 
to Hm, in which less than one-half the operators are of order 2 it must have 
such a co-set with respect to every possible Hm. Hence the possession of such 
a co-set is an invariant property of the group. When r= 2 this is clearly 
not the case, but when this condition is satisfied m is an invariant, being equal 
to unity, and H, is one of two groups. One of these must be of index 2 while 
for the other r may exceed 2. 

If both of these subgroups are of index 2 exactly one-fourth of the opera- 
tors of G have an order greater than 2 and G is known to be the direct product 
of the octic group and an abelian group of order 2% and of type (1,1,1,°-- ). 
When one of them has an index greater than 2 then each of the co-sets which 
arise with respect to this subgroup has exactly half of its operators of order 2 
while the remaining operators of this co-set are of order 4 and have a common 
square. The co-set with respect to the other of these two subgroups contains 
more operators of order 4 than of order 2, the number of the operators of order 
2 being half the order of the smaller subgroup which may be used for H,. 
The least order which this subgroup can have when m’ is odd is $(m’ —1). 
When m’ is even it is $(m’—2) since the commutator subgroup appears in 
the central of G. 

In this system of groups for which r may be equal to 2 it is evident that 
there is an infinite number of groups in which the number of operators is 
equal to the order of the group multiplied by a number of the form Ala 
where n is an arbitrary positive integer. In fact, if one such group is found 
we can find an infinitude of others by forming the direct product of this 
group and abelian groups of order 2¢ and of type (1,1,1,--- ), where a 
represents a positive integer but is not otherwise limited. 
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§ 3. Exactly one-half of the operators of each co-set are of order 2 and m > 1 


When m > 1 and exactly one-half of the operators of each co-set are of 
order 2 it is easy to see that it is not possible that the subgroup of index 2 
under H» composed of all its operators which are commutative with each of 
the operators of one co-set is the same for all the different co-sets. Let ¢ be 
any operator of order greater than 2 contained in G and let s, be any operator 
of Hm which is commutative with ¢ but not contained in the central of G. Let 
H,, as before, represent the subgroup composed of all the operators of G 
which are commutative with s,. 

The product obtained by multiplying any operator of order 2 contained 
in H, into an arbitrary operator whose order exceeds 2 in G— H, is trans- 
formed into its inverse by s, whenever the order of this product exceeds 2. 
In this case the two factors of this product must be commutative. As every 
operator whose order exceeds 2 in G— H, is transformed either into itself 
or into its inverse by every operator of order 2 in H, it results that all the 
operators of H, transform every operator of G— H, whose order exceeds 2 
either into itself or into its inverse. The squares of each of the operators of 
H, is therefore contained in the central of G. In particular, ¢? is in this 
central. As ¢ was any operator of order greater than 2 contained in @ it 
results that the squares of all the operators of G are found in the central of G 
and hence G involves no operator whose order exceeds 4. 

Since the quotient group of G with respect to its central involves only 
operators of order 2 besides the identity this quotient group must be abelian. 
In particular, both H, and Hm must be invariant subgroups of G. Every 
non-invariant operator of order 2 contained in G@ has just two conjugates. 
For, if such an operator s, had more than two conjugates, it may be supposed 
to appear in Hm. It may be assumed that s, is transformed into itself multi- 
plied by the three distinct operators so, s’9, S) 8’) of order 2 contained in the 
central of G. Let t), t’o, tot’, be three operators of G which effect these trans- 
formations respectively. If ¢, and ¢’, were not commutative with the opera- 
tors of the same subgroup of index 2 under H» there would be a co-set cor- 
responding to Hm which would have less than half of its operators of order 2. 

Hence it may be assumed that ¢, and ¢’, are commutative with the opera- 
tors of the same subgroup of index 2 under G. There must be an operator ¢’’, 
in G which is not commutative with all the operators of this subgroup. The 
operators of Hm which are commutative with both of the operators ¢, and ¢/’, 
constitute a subgroup of index 4 under H» and include the central of G. 
Hence ¢, and ¢’’, are non-commutative with the operators of a common co-set 
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with respect to this subgroup and it may be assumed without loss of generality 
that they transform the operators of this co-set into themselves multiplied by 
two distinct operators since ¢, may be replaced by ¢’). Hence we are again 
led to a co-set with respect to Hm which involves more operators whose orders 
exceed 2 than operators of order 2. It has therefore been proved that when 
m > 1 for every possible choice of s,, and one-half of the operators of each 
co-set with respect to Hm are of order 2 then each non-invariant operator of 
order 2 contained in G has exactly two conjugates under G. 

From the fact that each non-invariant operator of order 2 contained in G 
has exactly two conjugates under G it is easy to determine additional funda- 
mental properties of G. In particular, it is possible to select a set of inde- 
pendent generators of G in such a way that all of these generating operators 
are of order 2 and that each one is commutative with all of the others with 
the exception of at most one of them. In fact, we may select s, as the first 
independent generator of G and then let s, be any operator of order two in 
G—H,. The operators common to the two subgroups composed of all the 
operators of G which are commutative with s, and s, respectively constitute a 
subgroup of index 4 under G which includes the central of G. The order of 
the product of s, and s, is 4. . 

If this subgroup of index 4 is abelian it must be of type (1,1,1,--° ) 
and H, must be abelian and of the same type. In this case m would be equal 
to 1 and hence we may assume that the subgroup in question is non-abelian. 
More than one-half of its operators must be of order 2 since more than one- 
half of the operators of G are of order 2. It is therefore possible to find two 
additional non-commutative generators of order 2, s, and s,, such that each 
of these generators is commutative with both of the operators s,, s, and that 
(Ss S,)? = (8, 8.)%. If the last condition were not satisfied s,,s, would be an 
operator of order 2 having more than two conjugates under G. 

When the group composed of all the operators of G which are commu- 
tative with each of the four operators s,, s., 83, $, is abelian and of order greater 
than 2 it is of type (1,1,1,--~ ) and G is the direct product of a subgroup 
of index 2 under this abelian group and the group of order 32 generated by 
the given four operators. On the other hand, when the operators cf @ which 
are commutative with each of the four operators s,, s,, 83, $, are not all com- 
rutative with each other then the subgroup generated by these operatzrs must 
involve two non-commutative operators of order 2, s, and s,, such that the 
square of their product is equal to (s, s,)?. 

This process can evidently be continued until we obtain either a set of 
generators of G or an invariant subgroup of G of order 27+ such that @ is 
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the direct product of this subgroup and an abelian group of type (1,1,1,--- ). 
Hence we have established the following theorem: When m exceeds unity for 
every possible choice of Hm and one-half of the operators in each co-set of G 
with respect to Hm are of order 2 then G is generated by B operators of order 
2 such that each of them ts commutative with all of the others except at most 
one of them and the square of the product of every pair of non-commutative 
operators of this set is the commutator of order 2 contained in G. 


§ 4. Hzactly one-half of the operators of each co-set are of order 2 and m=1 


Two cases require consideration according as m—1 for every possible 
choice of operators from G for s, or as m = 1 for some one but not for every 
possible such choice. We shall first restrict our attention to the former of 
these cases and hence it results that the subgroup of index 2 under H, com- 
posed of all the operators of H, which are commutative with the operators of 
one of the r—1> 1 co-sets constitute the central of G. 

Let ¢, and ¢, be two operators of order 2 found in any two distinct co-sets 
with respect to H,. If ¢, and t, were commutative the central of G together 
with ¢, and ¢, would generate an abelian group of type (1,1,1,:-~- ) and of 
a larger order than the order of H,. By using ¢, instead of s, we would 
therefore find a new group for H, which would give rise to a co-set in which 
less than one-half of the operators are of order 2. As this case has been 
considered it may be assumed that ¢, is not commutative with any operators 
of order 2 except those found in the central of G or in the co-set to which ¢, 
belongs. 

Since ¢, transforms into its inverse each of the products obtained by 
multiplying ¢, successively into all the operators of order 2 contained in G, 
and as these products include all the operators of G whose orders exceed 2 it 
results that every operator G whose order exceeds 2 is transformed only into 
itself or into its inverse under G and that it must be transformed into itself 
by all the operators of G whose order exceeds 2 since it is transformed into its 
inverse by a number of operators of order 2 equal to one-half of the order of G. 
Hence it results that the operators whose orders exceed 2 in @ generate an 
abelian group and that G is either the dihedral or the extended dihedral group. 

It remains to consider the case when r > 2 and H, involves an operator 
s, which transforms every operator whose order exceeds 2 in & into its inverse 
and when, moreover, G contains a non-invariant operator s, which is commu- 
tative with operators whose orders exceed 2 in G. If the central of @ is 
composed of half of the operators of H, then s, must be contained in an 
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abelian group of type (1,1,1,- - - ) whose order exceeds that of H,. As one 
of the co-sets of G with respect to this abelian group would contain more 
operators whose orders exceed 2 than operators of order 2 it may be assumed 
that the order of the central of G is less than one-half of the order of H,. It 
therefore results that the square of every operator of G is found in the central 
of G and hence H, is an invariant subgroup of G. 

We may form Hm, m > 1, by, starting with s, instead of with s,. Let 
H’, be the subgroup of G composell of all its operators which are commuta- 
tive with s,. The product of an operator t, of order 4 in H’, and an operator 
t, of order 4 in G— H’, is transformed into its inverse by each of the opera- 
tors s, and s, whenever the order of this product exceeds 2. In this case, it 
results therefore from a transformation by s, that ¢,-* = and hence 
t, and ¢, are commutative. On the other hand, it results from a transforma- 
tion by s, that ¢,¢,-1 = ¢,1t,1. As the latter is impossible it follows that the 
product of an operator of order 4 in H’, and an operator of this order in 
G — H’, must be of order 2. 

As t, is commutative with one-half of the operators of Hm and is trans- 
formed into its inverse by the other half it results that it is commutative with 
all the operators of order 2 in H’, — Hm. Since f, is any operator of order 4 
in G— H’, all of these operators are commutative with the operators of the 
same subgroup of index 2 under Hm. This is contrary to the hypothesis that 
S. is commutative with ¢, but not with ¢,. Hence it results that if exactly 
half the operators in each co-set with respect to every possible Hm are of order 
2 and m=—1 for one choice of s, it must be equal to unity for every such 
possible chotce. Hence the central of G is of index 2 under H,. 


§ 5. Conclusion 


The most familiar system of groups in which at least one-half of the 
operators are of order 2 are the dihedral and the generalized dihedral groups 
and when the order of a group is divisible by an odd prime number there are 
no other groups which have the property that at least one-half of their opera- 
tors are of order 2. When the order of G is 2” there are other infinite sys- 
tems of groups’such that the number of operators of order 2 in the group 
exceeds half the order of the group but in all these groups the order of each 
operator is a divisor of 4. The dihedral and the generalized dihedral groups 
may be characterized by the facts that at least one-half of their operators are 
of order 2, that all the operators of such a group which are commutative with 
an arbitrary non-invariant operator of order 2 constitute an abelian group of 
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type (1,1,1,--- ) and that they do not contain an abelian group of this 
type which is of even order and of index 2 unless the central is of index 4 
or of index 1. 

One of the simplest systems of groups which are neither dihedral nor 
extended dihedral but have the property that more than one-half of their 
operators are of order 2 is obtained by extending the abelian group of order 
2-1 and of type (1,1,1,-°.-- ) by means of an operator of order 2 which 
transforms this abelian group into itself but is commutative with less than 
one-half of its operators. The order of the central of such a group cannot be 
less than 4(m’—1) since the multiplying group must appear in the central. 
Hence these groups are possible when and only when m’ > 4 and the number 
of such groups for a given value of m’ is $(m’—1) —1 when mm’ is odd and 
4(m’ — 2) —1 when m’ is even. It may be noted that for each of these 
groups m = 1, but that r has two possible values. 

As a second system composed of groups which are neither dihedral nor 
extended dihedral but have the property that more than one-half of their 
operators are of order 2 we may consider those for which one co-set with 
respect to every possible H» has fewer operators of order 2 than of larger 
order. The smallest group of this system is of order 64 and every other group 
of the system is the direct product of this group of order 64 and an abelian 
group of order 2° and of type (1,1,1,:-.- ). Exactly 7/16 of the operators 
of each of these groups are of order 4. 

Finally, there is a third system of groups which are neither dihedral nor 
extended dihedral but contain more operators of order 2 than operators of 
larger order. ‘This system is characterized by the facts that half of the opera- 
tors in each co-set with respect to Hm are of order 2 andm>1. The smallest 
order of such a group is 2° since if it is not dihedral or generalized dihedral 
it must contain at least two pairs of non-commutative generators when it 13 
generated by a set of B = 4 operators of order 2 such that each operator of 
this set is commutative with all except at most one other operator of the set. 
Every non-invariant operator of such a group has two conjugates under the 
group. By adding to the three systems just described the system composed 
of the dihedral and the extended dihedral groups we obtain the four possible 
systems of groups of order 2” which involve all the groups whose order is of 
this form and which have the property that more than one-half of their opera- 
tors are of order 2. 

To exhibit an important analogy between the second and the third system 
described above it may be desirable to note that the second system may be 
characterized by the fact that each of its groups may be generated by a set of 
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8B >S4 operators of order 2 such that each of these operators is commutative 
with all of the others except at most one of them and that there are exactly 
two pairs of non-commutative operators in the set, whose products are of order 
4 and have different squares. This system can also be defined by the fact that 
each of its groups contains the direct product of two octic groups, and when its 
order exceeds 64 it is the direct product of this group of order 64 and an 
abelian group of type (1,1,1,:-- ). 


— 
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The Green’s Function for a Plane Contour. 


By Hopart Dickinson FRARY.- 


A Green’s function is a function of a boundary and two points. The 
points, of which one is considered fixed and the other variable, both lie 
inside (or both outside) the boundary. In two dimensions the boundary 
is a contour. A large part of mathematical physics is devoted to the solution 
of boundary problems in two dimensions. All such problems can be solved 
immediately if the Green’s function for the given boundary or contour is 
known. 

The character of the Green’s function depends on three things,—(1) 
the differential equation of which it is a solution, (2) the contour or 
boundary to which it applies, and (3) the boundary conditions. The class 
of Green’s functions to which we shall confine ourselves in this article is 
the one of most frequent occurrence in mathematical physics, namely, the 
Green’s function for Laplace’s differential equation, with what Hilbert calls 
boundary condition (I) ; that is, G0 on the boundary. Having found this 
particular solution G, we can solve the problem of Dirichet in the plane; 
that is, we can determine a solution w satisfying any given condition of the 
form u=U(6) on the boundary. This we can do by making use of the 


formula 


which in polar co-ordinates becomes 


u(r,6) U (6) + | 


There are two well recognized methods of obtaining the Green’s function 
in two dimensions, both of which, however, are very limited in their applica- 
tions. These methods are the following: 

(1) The method of images. This method is easily applicable to the 
circle, semi-circle, infinite strip, half plane, equilateral triangle, 30 degree 
right-triangle, and the rectangle; but for other contours it would involve 
complicated Riemann surfaces and hyper-elliptic functions. Apparently no 
one has attempted to apply the method to other polygons. 

(2) The method of Schwarz. Schwarz obtained a general formula 
w = f(z), where f(z) is in the form of a definite integral, by means of which 
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any given polygon can be mapped conformally on to the unit circle. Knowing 
this mapping function f(z), we can obtain the Green’s function for the poly- 
gon immediately from the formula: 


= — Re log f(z) (Re = real part of) 


By this method we are able to get the Green’s function for the equilateral 
triangle, the 45 degree right-triangle, an infinite strip, a regular five-pointed 
star, and a rectangle. Beyond this it is difficult to go. Other polygons lead 
to Abelian integrals and hyper-elliptic functions. Moreover, except for regu- 
lar polygons, it seems to be impossible to determine the unknown constants 
which appear under the integral sign in the function f(z). For regular 
polygons, however, we could evaluate Schwarz’ integral, 


dz 
w= 9 
(1—2")" 


in the form of an infinite series, and take the real part of the logarithm. In 
other words, the series obtained by expanding 


é 08 (1—2")” 
gives us the value of the Green’s function for a regular polygon of n sides in a 
circle of radius 


dz 


where the point of discontinuity is at the center of the polygon. This work 
of Schwarz has been considerably discussed and amplified by other writers, 
particularly in a recent book by Study.* However, very little of fundamental 
importance has been added to Schwarz’ original memoirs. 

The method of arriving at the Green’s function which we shall employ 
in this article is new, although in principle similar to the methods used by 
Fourier and Neuman in attacking various physical problems. We select an 
infinite set of linearly independent functions which are solutions of the given 
differential equation, and seek to expand the Green’s function, or rather, the 
quantity G—y, where y is the principal solution (log 1/r in this case), in a 
series of these functions, determining the coefficients of the series so that the 
boundary conditions are satisfied. The determination of these coefficients 
depends upon the solution of an infinite set of equations in an infinity of 
unknowns. In the problems solved by Fourier the solution of the infinite 


* Vorlesungen iiber Ausgewihlte Gegenstinde der Geometrie, 2tes Heft. 
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set of equations was simple. In the problem before us the solution is much 
more involved, though entirely manageable in a large number of cases. 

This method of procedure is much more widely applicable than either 
of the two methods described above. These two methods give workable results 
for, at most, the regular polygons, two other polygons, the five-pointed star, 
and the rectangle. Our method is applicable, not only to all regular polygons, 
but to a large class of irregular polygons as well, and to many other contours. 
It gives the Green’s function in the form, log 1/r plus an expansion in cir- 
cular harmonics, r" cos n@ and r” sin n6. The coefficients of this expansion 
are each given in the form of an infinite series of constants and can therefore 
be determined to any desired degree of accuracy by taking a sufficient number 
of terms. We are convinced that the method is actually workable in a practical 
sense in a large number of cases. 

We have confined our attention to a single class of Green’s functions, but 
the same procedure might be employed to get the Green’s function for any 
differential equation with any given boundary conditions.* It might also be 
extended to problems in three dimensions. 


Formulation of the Problem. 


By a plane contour we mean any closed plane curve which is regular in 
the sense of Osgood; that is, is composed of a finite number of analytic arcs 
or straight segments. An important case is that of the polygon. 

The Green’s function for any such contour is given by the following 
formula in polar co-ordinates: 

G(r, 0) = —logr + u(r, 6), (1) 


where the function wu is harmonic, that is, satisfies Laplace’s equation, 
low 


and has no singularities in the region. 

Consider this function, u(r,6). We know that it is everywhere harmonic 
in the area S enclosed by the contour C and has no singularities in 8. More- 
over, on the boundary, for which r— and u—U (6), we have G = 0 and 


U (R, 6) =log R(6). (2) 


Draw any circle, radius r, (Fig. 1), about 0 so that it lies wholly in 9. By a 


*See Hadamard’s Columbia lectures, pp. 47-52. 
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well-known theorem we can expand the harmonic function u uniquely in the 


series, 


(An cos nO + By sin 6), (3) 
1 


Where A, and By are the Fourier constants for the area enclosed by the 
circle: OS 6S 22,0Sr=r,. The functions r” cos n6 and r® sin n6 form 
a complete set of orthogonal functions for the circular boundary. Let us write 
equation (3) for convenience in the form 


8), (4) 


where ¢, = 1, ¢, = rcos6, 


Fre. 1. Fie. 2. 


We can continue the function w analytically across the boundary C (Fig. 
1) by means of Schwarz’ principle of symmetry. The extended function thus 
obtained takes values equal but opposite in sign at every pair of points which 
are symmetrical * with respect to an arc or segment of the contour. At the 
symmetrical points or images 0’ of the point 0, the function becomes infinite 
as log r, but at all other symmetrical points it is harmonic. Hence if we 
construct the images ¢ of the point 0 in the several arcs or segments of the 
contour, we know that the function u exists and is harmonic at least in a circle 
extending to the nearest image. Fig. 2 is drawn for a special case where the 
contour is a polygon. 

Now, the question arises: Is the function w represented by the same 


* For definition see Osgood, p. 671. 
7 For the method see Lery’s thesis on La Fonction de Green pour une Contour 
Algébrique, p. 49. 
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expansion in the whole circle? To answer this question we cite the following 
theorem due to Bécher.* 


Theorem: If uw ts a function harmonic throughout the neighborhood of 
the point xy, and if, when this function is continued analytically, the distance 
from XYo to the nearest singular point of u, which lies in the same sheet of 
the Riemann surface generated by the analytic continuation in which x,y, lies, 
is K, then the development + in polar co-ordinates (the point roy, being taken 
as pole), 

u=A,/2 + (Anr” cos n6 + sin n6), 


converges and represents u throughout the interior of the cwrcle of radius 
r = K and does not converge throughout any continuum which does not lie 
in this circle. 

If we restrict the problem, then, to contours of such a nature that a point 
0 can be found whose images all lie outside of a circle about 0 enclosing the 
contour (this we can make the unit circle), we have from the above theorem, 
for all contours of this type and all such points 0 within each contour, the 
Green’s function for the point 0 and the given contour in the form 


G=—logr+ A,/2 +3 cos n6 + Barr” sin 
1 


= — log + 6). (5) 


This expansion holds for all points inside the polygon and on the boundary. 
Moreover it is unique and we must be able to determine the constants An and 
B,, or Cm from the boundary conditions. 

For the boundary we have r= R(6) and G =), that is, 


A, /2 + 3[AnR"(8) cos n6 -+ BuR"(8) sin n6] = log R(8) 
1 
co 
or > log R(6), (6) 
0 


where ® = 1; ®yn_, = R" cos nb; = sin nd. 


Now we can apply a method of analytic continuation due to Borel, to 
obtain a new expression for the function w which will hold for a much more 
extended region, as follows: Equation (3) may be written 


co 
A,/2 Re — iBn) 
1 


* Trans. Am. Math. Soc., Vol. 10, p. 278. 
+ Bocher gives here also an equivalent in terms of polynomials in # and y. 
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where 2 is a complex variable. Since uw is harmonic the expression 
oo 
>(An—1B,) 2" is analytic in a circle extending to the nearest image 0’ 
1 
of 0. Consider for example the polygon * ABCDE (Fig. 3) with the point 


of discontinuity at a point 0 inside it. By means of Borel’s integral we can 


oo 
continue the function = Anz" analytically to the region enclosed by the polygon 


A’B’C’D’E’ whose sides are lines through the several images 0’ and perpen- 
dicular respectively to the several lines 00’.. The expression for the function 
thus extended is 


0 


Applying this to the problem in hand, we have, in the real plane, 


oO 
A,/2 + Re (2t)" 
0 1 


co fe jn 
= A,/2+ fey (Anr® cos n6 + Byr” sin n6) dt. 
n! 
0 1 


This function is harmonic inside the polygon A’B’C’D’E’ whose sides are 
respectively parallel to the sides of ABCDEH. Hence this gives us a formula, 


fo fo @) jn 
G=—logr+A,/2+ Set} (Ane cos nd + Bur® sin nb) — dt, (7) 
0 1 


*Some or all of the sides of the polygon could be replaced by analytic arcs. 
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for the Green’s function for any convex polygon whatever, with the point of 
discontinuity at any point whatever inside the polygon. Evidently this 
method can be applied, not only to polygons, but to any convex contour what- 
ever and the point of discontinuity may be taken at any point whatever inside 
the contour. 

The whole problem is then reduced to the method of determining the 
infinite set of constants A» and Bn, or Cm. This problem may be solved in 
two ways: (1) by integral equations and (2) by an infinite set of algebraic 
equations. The first method gives only a formal solution; the second is a 
practical method as we shall show. 


Formal Solution by Integral Equations. j 
The c’s are Fourier constants. Hence, by a well-known theorem of 


co 
Fischer and Riesz, the sum of their squares 3 c¢?m, converges. Moreover, if 
0 


we select any set of orthogonal functions whatever, linear, planar or solid, 
for example hm(x), there exists one and only one function H(x) for which Cm 
are the Fourier constants with respect to this orthogonal set. Whether the 
series thus obtained converges or actually represents the function at all is 
immaterial; for in any case the formulas for the coefficients hold, namely: 


Cm = SO) hm(x) dx, (8) 


where the definite integral is either single, double or triple, according as the 
orthogonal functions chosen are linear, planar or solid. Now substitute these 
constants back in equation (6) and we have, provided the series in the bracket 
converges, 


f [y Rim (x) m®m (4) | dy = log R(6), (9) 


an integral equation of the first kind with unsymmetric kernel K(x, 6) = 
oo 
him (x) (8). This integral equation has been studied by Picard, 


Schmidt, Bateman and others, and existence theorems have been worked out; 
but no general method of solution has yet appeared in the mathematical 
literature. The following special form of the equation, however, is of con- 
siderable interest : 

If we take for the set of orthogonal functions hm(x) the linear set cos ny, 


2 
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sin nx, the function H(x) is nothing else than the value of w on the unit 
circle, 


A(x) =u(1, x); 
and the kernel reduces to the well-known expansion for 


1— 
1— 2 Rcos(x— 6) + R? R<1 


Thus the integral equation (9) now becomes 


23 1— R?(6) 
1/2 f, x) cos(x — 6) 


where the left hand member is Poisson’s integral. Hence we may state the 
solution of the problem before us as follows: 
Consider a closed contour of such a nature that a circle, with center at 
a point 0, can be drawn completely enclosing it, but not enclosing any of the 
images of 0 with respect to the contour or arcs composing the contour. Make 
this the unit circle. Then the Green’s function for this contour is obtained 
by using Poisson’s integral backwards, as it were, to determine the value on 
the unit circle, namely H(x) = u(1, x), of a harmonic function u(r, 6) which 
reduces to log R(@) for points of the given contour. Having found H we 
determine the Fourier constants or coefficients corresponding to it and substi- 
tute them in formula (5) to get the Green’s function. 
On account of the complicated character of the kernel in equation (10) 
the solution of this equation seems to be impracticable. We shall therefore 
take up in the following pages an entirely different method for determining 


the constants for equation (5.) 


—=logR(8), (10) 


Practical Solution by an Infinite Set of Equations. 


Going back to equation (6), it is possible to reduce the problem of deter- 
mining the unknown constants cm to the solution of an infinite number of 
equations for this infinite set of unknowns. One way to do this would be to 
expand the various powers of R(6) as well as the function log R(6) in Fourier 
series and compare the coefficients of cos n@ and sin n@ on each side of the 
equation. This leads immediately to an infinite set of equations for the 
constants Cm, but these equations are very awkward to handle. 

In order to get a more easily manageable set of equations we will con- 
struct a normalized, orthogonal set of functions from the set 


R(6) cos 6, ®, = R(6) sin 6, ®, = R?(6) cos 26, 


| 
! 
! 
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functions which will be orthogonal for the variable 6 between 0 and 2z. 
These functions we can construct by the method of Goursat, provided the 
functions are linearly independent. We shall first investigate this matter of 


linear independence. 
An infinite set of functions ®,(6) are said to be linearly independent if 


no finite number of them are linearly dependent; that is, if no equation of 
the form 


+ Pr, + + + CnPr 0, (11) 


where the c’s are constants not all zero, exists between a finite number of them. 
For the set of functions we are considering there certainly do exist such rela- 
tions for some contours r= R(6). For example, if sec we have 


&,— o,=0; 26, —O,—0. 


However, for no elosed finite contour can such a relation exist. For, suppose 
we consider a contour r= (6) such that a linear relation of the form (11) 
exists. Then the left hand member of this equation is the value V(R,6) of a 
harmonic function V(r,6) on a finite closed contour; hence by Gauss’s mean 
value theorem, 

V(r, 6) =0 


for all points inside the contour, which is only possible if all the coefficients 
vanish. Hence the set of functions 


&, = 1, 6, Roos 6, &, = Rsin 6, ©, = R? cos 26, 


and in fact any infinite set &,(@) obtained in the same way from a linearly 
independent set of two-dimensional harmonics ¢n(r, 6), are linearly inde- 
pendent for all closed contours, in particular for polygons. 

We can now apply Goursat’s method and form a set, %, ¥,...... , 
of normal, orthogonal functions, as follows: 


Vy = Ago Py 

VY, = + 

VY. = Ago Py + + 

VY; = Ago Py + + + (12) 


Vn = AmoPy + Ami P, + + + + AmmPm 
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Solving for the ®’s, we have 


B= 1 
= R cos 6 + b,,¥, 
b, = sin 6 = + + 


Don, = cos NO = dona, Vo bona, bona 
+ Don-1, on-1Von_1 


(13) 


+ Don, on-1 + Don, onVon 


where 
Ams, m mer 0 
Ams2, m Ams2, ms2 0 

1 
1)* 
Om, m Om41, m4. °°" Ann 

An, m On, Un, 

Now expand log R(@) in terms of %, thus: 


co 
log R(6) 


where the coefficients are given by 


(14) 


27 
Ar= log R(6) a6. (15) 
0 
Substituting in equation (6), we have 
co 
=0 
where Co, Ci) Coy °° °° * are the unknown constants required. Comparing 
coefficients of the ¥’s on each side of the equation, we have 
Boo Co B19 Cr + Co + Cg 
Day Cy Og =A, 


? 


an infinite set of equations for the c’s. Solving this set of equations formally 
for each c by eliminating the other c’s one at a time, we have an expression 


| 
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for each c in the form of an infinite series involving the b’s and the A’s. If, 
then, we substitute for the b’s from equation (14), we arrive at the simple 


result 


fo 
Cm = 
i=m 


Substituting this solution in equation (6), we have 


+ (3 as, on, AiR" cos nb + di, on AiR" sin nO) = Ai. 

4=0 n=1 4=2n-1 4=2n i=0 
It is sufficient to multiply equations (12) by Ay, 4,, 4. °°° respectively, 
and add, to verify formally the correctness of this equation. Now since we 
know the c’s exist and are unique, if we can prove that equations (18) give 
an actual solution of (17), this is the only solution. To prove this we shall 
make use of the following theorem due to Anna J. Pell: * 

Theorem: If the sequences {Ax} and {ui} are such that the sequence 
{bix/Ac} is of finite norm for every &, and the matrix (Ajax: /yi) is limited, 
then for every sequence {Az} such that {,Ax} is of finite norm, the system 
of equations 


| (k= 0,1, 2, ) 
4=k 


has a solution c, such that {Ax¢c.} is of finite norm, and the solution is given by 


co 
Cy = 0, 1, 3, 3,° ) 


Making x = »; = 1, we can show that the three conditions required by 
the above theorem are satisfied, as follows: 

(1) The first condition to be satisfied is that the sequence {b4x} shall be 
of finite norm for every &; that is, the coefficients in equations (13), taken by 
columns, must give sequences of finite norm. This we can prove as follows: 
Since the &’s are normal, orthogonal functions we have, for every n, 


* Ann. of Math., Vol. 28, p. 35. We have changed Mrs. Pell’s notation to conform 
to our own. 
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27 
27 


Hence for every k, we have 


27 27 


2n-1, k 
2T 27 
R" sin n6 dd < Se d6 = Rn 
so that we have, for the norm of the sequence {Dix}, 
co co 
k=0, 1,2 (20) 


But the maximum value of (6) is always less than unity. Hence the right 
hand member of equation (20) is an absolutely convergent series. Thus we 
have proved that {biz} is of finite norm for every k. 

(2) The second condition is that the matrix (axi) of the coefficients in 
equations (18) shall be limited. This we can show to be true as follows: 
Consider first the matrix (bi) of equations (13). From (19) we have 


4 co 
302. Rh" <M, a constant; 
i=o Man 


that is, the b’s are of finite norm and hence * the matrix (bi) is limited. 
Now the matrix (ai) of equations (12) is the unique reciprocal of the matrix 
(bi), and hence is itself limited.+ Therefore its conjugate, the matrix (ax:) 
of equations (17) is limited. 

(3) The third condition is evidently satisfied; for the sequence {Ai} of 
the right hand members of equation (17) is of finite norm since the A’s are 
the Fourier constants for the expansion of log R(6) in terms of the W’s. 

Now by Mrs. Pell’s theorem it follows that equations (18) give an actual 
solution of equation (17), and hence give us the required expressions for the 
constants c. 

Formula (5) for the Green’s function may now be written 


G =— logr + + + Awe 
where 


Yo = Ao 


* Hellinger and Toeplitz: Math. Ann., Vol. 69, p. 307. 
Hellinger and Toeplitz: Math. Ann., 1. p. 311. 
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Ws = Ago + COS O Sin 6 + cos 26 


Now substituting for the coefficients A; from the equation (15), we have 
co 27 
G—=—logr 3 yi(r, 6) flog dd, (21) 
i=0 0 


which, since the matrix (ai) is limited, can be written 


=—logr+ 3 6) dé. (22) 


Thus we have an expression for the Green’s function in terms of a line integral 
around the boundary. When r= R(6), yi(r, 6) = (6), and the last mem- 
ber of the equation (21) becomes the expansion for log F in terms of the W’s, 
so that G vanishes identically. 


Application to a Square Contour. 


As a simple application of the above, let us find the Green’s function for 
a square with discontinuity at the center. If we make the apothem of the 
square unity and measure angles from it, we have for the function R 


R(6) = sec 6 —r/4= = 
R(6) = csc 6 =3n/4 (23) 
=— esc 6 = /4. 


Beginning with the first function 6,—1 we shall first orthogonalize all the 
other functions to it by determining OC _, and C to satisfy the following 
equations : 
27 
(Don, — B dd = 0 (24) 
27 
Soo (Gan — 0. 


From the first of these equations we have Cyn, 0 for n odd and n/2 odd; 
that is, 

For n/2 even 


Con = 2/2 f sec"6 cos n6 dé. 
/4 
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From the second of equations (22) we have Cyn—0 for every n. Thus we 
obtain a new set of functions ®’, all of which are orthogonal to ®,, as follows: 


= cos 6 sin 6 

= RF? cos 26 = sin 20 
= cos 36 sin 30 

&/ cos 46 — .48826 = sin 46 

& = cos 56 = sin 56 
cos 86 — 1.11259 & — R sin 86 

&/ — R" cos 120 — 3.01549 & — sin 120 

Now we take ® and orthogonalize /, ®{, ®,:-:-:-: to it in the same way 


as before. This process is continued until a sufficient number of orthogonal 
functions has been obtained. Finally these orthogonal functions must be 


normalized by dividing each function ® by [ fo ao |. Denote by %, ¥,, 
0 


Resi , as in the preceding section, the normal, orthogonal functions 


obtained in this way. _ 
Next we determine the A’s by means of the formula 


27 
n= f log dé. 
0 


Finally, substituting in equations (18), we have the required constants for 
the square, namely, 


C= A; = — .23653 —.05516 +.00883 +.00474 +--+ == 
C, =C,=0 

—,11296 —.02100 —.00834 = —14-+ 
Cy = Co = = Cy, = = O13 = = O 

= —.00794 —.00013 + ---—=—.01::: 
Cyz = Cyg = = Coy = Coy = Cog 0 

Cos +.00153 + +++ = 


Putting these constants in equation (5), we have the following formula for 


| 
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the Green’s function for a square with the point of discontinuity at the center 
and diagonal of length 2: 


= — log r —.28 —.14 r* cos 46 —.01 r* cos 86 +-.002 r** cos 126 +--:-- 
= Re (— log —.28 —.14 —.01 28 +.002 4+ ) (25) 


approximately. 

Thus we see that the expansion for wu, that is, G+ logr is the real part 
of an expansion in fourth powers of z. This may be verified by examining 
the expressions in elliptic functions, obtained by the method of images, for 
the square with the point of discontinuity and the origin of co-ordinates at the 
center. The character of the expansion for the Green’s function for the square 
suggests immediately that for a regular polygon of p sides (p any positive in- 
teger) with discontinuity at the center, w will be the real part of an expansion 
in pth powers. It is evident at once that this is true when p becomes infinite, 
for in that case we have a circle for which wu is constant. This may be shown 
to be true in general by making use of Schwarz’ mapping function for a 
regular polygon of p sides. 


“4 
| 


| 


On the Solution of Certain Types of Linear Differential 
Equations in Infinitely Many Variables. 


By WEBSTER G. SIMON. 


The main purpose of this paper is to prove the existence of certain types 
of solutions of particular kinds of linear differential equations with periodic 
coefficients in infinitely many variables. As a means to this end the existence 
of exponential solutions is established for certain types of linear differential 
equations in infinitely many variables with constant coefficients. 

The starting point is the existence theorem given by von Koch,* and a 
generalization of Poincaré’s theorem + concerning the development of the solu- 
tions of the differential equations as power series in a parameter » when the 
functions appearing in the differential equations are themselves power series 
in w. Then our work is very similar to that of the finite case, the finite 
determinants becoming infinite determinants which together with all their 
first minors converge absolutely. 

The type of determinant used in this paper is more general than the 
normal determinant, and is the following.§ The determinant 


1+ A125 Ay 39 
1 + Ago, 
Azo, 1 +- 


is such that there exist two sets of positive constants Sy 


which are of such a nature that | aij | < SiT;, and 3 S8iT; converges. This 


is a type of infinite determinant given by von Koch $ "which together with all 
its minors converges absolutely. 

Using the methods thus indicated, we find that many of the phenomena of 
the finite systems are carried over into the infinite systems of differential 
equations. | | 


* Von Koch, Ofversigt af Kongliga Vetenskaps Akademiens Férhandlingar, Vol. 56 
(1899), pp. 395-411. 

t Poincaré, “ Les Méthodes Nouvelles de la Mécanique Céleste,” Vol. I, chapter IT. 

§ Von Koch, Acta Mathematica, Vol. 24 (1901), pp. 89-122. Hereafter this paper 
will be denoted by K. 

||See Moulton and MacMillan, AMERICAN JOURNAL OF MaTHEMATICS, Vol. 33 
(1911), pp. 63-96. 
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I 


The differential equations which we shall consider are of the general type | 


00 
(1) = 64; (t) xj, 2, 


where 2 is the derivative of 2; with respect to the independent variable t. 
Concerning this system of differential equations we shall establish two 


existence theorems. 

Theorem I. Suppose the system (1) satisfies the following hypotheses: | 
(H,) The 6:; (t) are analytic functions of t when |t| < R. 
(H,) Positive constants, S:, +, exist such that | 6:;(t) | 
< SiT;,|t| < RB, and 8,T7,+8,T,+: converges. 

(H,) The (0) = Bi, where the Bi are constants such that B,T, + 
converges. 
Then there exists a unique system of functions 


co 
(2) Pi +2 += 1,2, 


which satisfy the system (1), and which converge for |t| < R. 


Although this theorem was established by von Koch,* it will be proved 
briefly here because the notation and the results are essential for the later 


parts of this paper. 
From (H,) we have 


Upon substituting the series (2) in (1) and equating coefficients, we get 
(1) (0) Q(0) 


co 
2p —3 + 


co 
(A 1) pat — > a where po Bi. 
j j $ 


. 


We see that the formal solution is unique. 


* Loc. cit. 
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To prove the convergence of (2), consider the system of differential 
equations 


(5) ji, t= 1,2, 


co 
Upon substituting = Bi + %y{/t in (5) and equating coefficients, we get 
j= 
Sid Ty, 
j=1 


co 
Qy2) 


co 
(A+ 1)yoarm == where = Bi. 
j=1 


On comparison of (6) and (4) we see that >| 14, ©. 
Hence & dominates x; for every 7. 
From (5) we have 


(7) 
On taking (0) = 0, we have 
(8) Bi=Sié, 1=1,2,°°* ©. 
Therefore each equation (5) reduces to 
(9) &=Cé§+K, 


where C= 

From the theory of a finite system of differential equations we know that (9) 
has a unique analytic solution for |¢| < R. On combination of this fact 
with (8), we see that the solution (2) of (1) converges when |¢| < R. 

Now we define as a fundamental set of solutions of (1), a set such that 
every solution of (1) can be expressed as linear homogeneous functions with 
constant coefficients of the elements of the set. Then if we denote by ¢i; the 
elements of the fundamental set, it follows that the determinant of the qi; 
converges and is not zero for all | ¢| < R, and conversely. Furthermore it 


Ze (t) dt 
has the form A= Ae? sit"? for all |t| <BR, where A, is the value of 


the determinant A when t= {¢,. If, for example, we take the solutions defined 
by =1 and =0, we see from (8) and (9) that 


7 
| pis | | —1) |, and oy —1) |, 


= 
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‘We see that the determinant of the ¢i; is absolutely convergent. 


Theorem II. Suppose that the system (1) satisfies the following hypo- 
theses: 
(H,) The 6:; are expansible as power series in a parameter p, which converge 
for all real |t| << Rif |p| <p. 
(H,) For p= 0, the 0:5 = aij, where the aij are constants. 
(H,) The 6:5 satisfy all the hypotheses of Theorem I uniformly with respect 
tomif|ul|<p. 

Then the solutions of (1) can be expressed as power series in p, which 
converge for | ~| Spo <p and for all real |t| < BR. 


co 
In Theorem I we showed that |t| <BR, |u| <p, 
j=1 
where np» = P® 


co 
In particular = ay The are power series in ». Hence 
j=1 


CO 
3 SB int which converge for <p. Then by 
a well-known theorem in the theory of double series, we infer that Bf) = 

3 By Ms |4|<p. In a similar fashion we proceed step by step, and show 


that the 8% are also power series in », which converge for |u| <p. So we 
have 
foe] 
(10) Li = dix (mu), += 1,2, 
k=1 
which is uniformly convergent for |»|Sp,<p, and where ¢ix(») = 


oy +40u+:°°-°. Then by a well-known theorem of Weierstrass, we 


know that we can rearrange the series (10), and write it in the form 
co 
=1 


II 


Now let us assume that the 6;; of Theorem I are constants; and let these 
constants be denoted by aij. We inquire whether the system (1) has a solu- 
tion of the form 


(11) 


co 
where by (H;) of Theorem I 3aijc; converges. Necessary and sufficient con- 


| 
| 
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ditions that such a solution exist are 


(12) 4 + —a) Co + + =O, 
In order to see that (12) can have a solution for which the c; are not all zero, 
we make the transformation a = — A’, assuming of course that a is not zero. 
Then the equations (12) become 
(1 + Cy Ady + 
(13) 4 + (1 + Adiga) Co + 


In order that the equations (13) have a solution in which the c; are not all 
zero, it is necessary and sufficient that A be a root of 


(1 ++ 


Equation (14) is called the characteristic equation. Since the ai; satisfy the 
hypotheses of Theorem I, it follows that A(A) is an integral function * of A. 

If A, is a root of (14), then there exists a finite number r such that 
A(A,) together with its minors of order 1,2, - - - , r—1, vanishes; but there 
is at least one minor of order r which is different from zero. Let the minor 
which is obtained by replacing the elements in the i,th row and the &,th 
column by unity and the remaining elements in that row and column by zero, 
and so on up to the 7-th row and the &,th column be denoted by 


tes 

Kez, 

and suppose that this minor is one which is different from zero. Then the 
system (13) admits the solution f 


(15) 


depending on r independent parameters &,,--*~- &, and admits no others. 
To sum up our results we have the following 


*See K, p. 104. 7 See K, p. 108. t See K, p. 109. 
P 


| 
| 
kee 
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Theorem III: If the 6:3 of (1) are constants and if A, is a root of 
A(A) =0 such that all minors of order r —1, but not all of order r, vanish 
for X=Ab, then there exist r independent solutions of (1) of the form (11). 


Now the characteristic equation, since it is an integral transcendental 
function, may have no finite roots, it may have a finite number of finite roots, 
or a denumerably infinite number of roots. In the first case the c; are all 
zero, in the second case, a finite number of solutions of the form (11) exists 
and in the last case an infinite number of the form (11) exists. This is the 
case we shall study, for this infinite set constitutes a fundamental set of solu- 
tions as the following discussion shows. 

For each A that is a root of A(A) —0 there is a solution of (1) of the 
form 


(16) = €450%7*, 4,9 1,2, °° ©. 


We know that the system (1) admits the fundamental set of solutions ¢4;, 
i,j=1, 2,° where $i:(0) —1, $1;(0) =0, 147. Consequently 


co 
k=1 
First let us assume that the a; are all distinct. Then it follows immediately 
that the system (16) constitutes a fundamental set. For consider 2i,, $j, 


$i» Pros 

(1) =A! Por» * * 
= Al Poo» 


Choose the initial conditions such that A‘ = 1, which is no actual restriction 
since we are merely determining the arbitrary constant which enters, and we 
can multiply our solution by a constant, not zero, when we are through if we 
choose. Thus dij, 1,2, 0, 2,3, 00, constitute a fun- 
damental set. Similarly we show that 2:,, 2:2, ¢i;, constitute a fundamental 
set, where 7 = 3,4,---+ o. Continuing in this manner and passing to the 
fimit we see that the wij, 1,7 oo, constitute a fundamental set. 
A result of this determination of the arbitrary constant is that ci; = 1. 

Next let a, a,, and aj;s4a,, 38,4,°°+ There are two cases 
that may arise here, viz., all of the first minors of A(A) vanish when a= a,, 
or not all vanish. Consider the former case. Then (15) assures us that we 
can solve for the ci; in terms of two of them. On choosing the notation so 


\ 
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that ( 4 is not zero, we get two independent solutions of the form . 


Liy = Dig = C1,0%', 1,2, - ©, 


which with 7 = 3,4, + + 0, constitute a fundamental set; for 


Pigs AM bi + AM 29 AM $12, Piss 
Loi» Loo» Pag = |AM¢,, + A bo, + 


and the A“, A may be taken to be zero, and A® =—A®=1. So these 
solutions, we see, constitute a fundamental set. 
Next let us assume that not all the first minors of A(A) vanish when 


a =a, = a;; and let us choose the notation so that the first minor G) is not 


zero. Then to get a solution associated with a,, we make the transformation 
+ and, if possible, deter- 
mine the b’s in such a manner that =ay. On these 
in (1) we get + + | + | 
=a[b,7, + ]. 
This equation must hold for all initial values of the 2’s. Therefore it me an 
identity in them; and the 6’s satisfy the following set of equations 


( (a,, —a)b, + =0, 
(17) ne + =O, 
In order that these equations have a solution for the 6’s not all zero, we get 
as a necessary and sufficient condition, just as in the case of c’s in (13), 
A(A) =0. This condition is satisfied, for we have assumed that when a —a,, 
the characteristic equation is satisfied. Furthermore we have chosen the nota- 


tion such that the minor ( 4 is not zero whena—=a,. Then we get for the b’s 
bj = mjb,, j=1,2,°-+ ©, 


1 
where mj; = () - For convenience we shall take b, —1. Then the equations 
1 


become 
y= = 141, 
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Next we make the transformation 
and determine, if possible, , such that = ay. + 
Then, as in the foregoing, we get the following, infinite set of equations for 
the determination of the d’s 


( (a, —a)d, + + + =1, 
(17) 0+ (a2 *=0, 
On setting a——1/A, we get an equivalent set of equations for the deter- 
mination of the d’s, the determinant of which is 
A, (A) = |0, 1 + — |= (1 + Aaj) — 


By making use of the equations (17), the reader will convince himself that 
A, (A) is identically equal to A(A), for it is obtained from A(A) first by inter- 
changing the rows and columns of A(A) and then multiplying the elements 
of certain rows by certain quantities and adding them to the corresponding 
elements of other rows. Furthermore there is a first minor of A,(A) which is 
not zero when a=a,. This minor involves the elements of the first row, 
since all the other first minors vanish. 

The determinant of the coefficients of the equations, omitting the first 
one, is zero for a= a,, since when a =a,, A(A) has a multiple root A,._ There- 
fore we can solve for the ratios of d,, d,,- ++. Then let us substitute the 
values of the d’s thus obtained in the first equation of (17’). First, we know 
that the sum will converge *; for the d’s gre proportional to the first minors 
of A,(A), and the sum of the products of the elements of a row and the cor- 
responding first minors converges. Secondly, since the d’s which we have 
determined carry an arbitrary factor, we can finally determine them so that 
the first equation is satisfied. 

Therefore our equations have been reduced to 


where the a® are the transformed ai; of the original equations. And on solv- 
ing the systems (1’) and (1) and putting in place of y, and y, the cor- 
responding value of xz, and z., we see that the solutions of (1) associated with 
a, are 

Li, = C4,0%*, == (Cin + )e%*,, 1,2, °° ©. 


* F. Riesz, “ Les systémes d’équations linéaires 4 une infinité d’inconnues,”’ p. 34. 
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Then by a precisely similar argument to the foregoing we can show that 
Viz, Vig, 1,2, ©, 3, 4,° °° 0, 


constitute a fundamental set. 
When for a given value of a the characteristic equation has a root of higher 
multiplicity, discussions similar to the foregoing must be made. 

We should expect to treat next the case that A(A)= 0 has a root of infinite 
multiplicity. But such a case cannot arise. For we have seen that A(A) is 
an integral transcendental function of A. ‘Therefore it can be expanded in a 
Taylor’s series in the neighborhood of any finite A. Suppose that A, were a 
finite root of infinite multiplicity. Then A(A) vanishes together with all of 
its derivatives when A=A,. Hence A(A) vanishes at every point in the 
neighborhood of A». Then it is identically zero. But this is not true, for 
when A= 0, A(A) =1. Hence A(A) has no root of infinite multiplicity. 

If a= 0, we can not make the transformation which carries the system 
(12) into the system (13), and the determinant of (12) diverges. Then we 
do not know whether the equations (12) have a solution for the c; not all 
zero; but each special case must be considered as it arises. 

To sum up our results we have 


Theorem IV. If and only if A(A) has an infinite number of roots ri, the 
system (1) has a fundamental set of solutions each of whose elements is of 
the form 


Lig = (t), 2, 


where a, = —1/d; and the Wi; are polynomials in t of degree at most (n—1), 
and n is the order of multiplicity of the root dj;. 


III 


Now we shall assume that the 6;; of Theorem I are periodic with the 
period 27. We have seen that the system (1) has, as a fundamental set of 
solutions, the set ¢i;(¢), where 


(18) ii (0) = 1, $i3(0) = 0, 1, = 1, 2, 
Let us make the transformation 
(19) erty, 


where a is an undetermined constant. Then the equations (1) become 


co 
(20) + ays = 3 yj, 1,2, 0. 


J=1 


r 
|| 
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Since the system ¢1; constitutes a fundamental set of solutions of (1), any 
solution of (20) can be written 


co 
(21) Ajpij(t), t= 1, 2,° 


We now inquire whether it is possible to determine the A; and a so that 
the y; defined by (21) shall be periodic with the period 27. From the form 
of (20) it is evident that necessary and sufficient conditions that the yi be 


periodic with the period 2m are 
(22) (2) — yi (0) = 0, t= 1,2, 


On imposing these conditions on (21), we get 


Then we make the transformation 
(24) —o/o —1, 


assuming of course that e°¢7=541. Then the equations (23) become 


oo 
(25) + As[1 + o (pis (2x) —1] = 0, 
j=l iA~j, i=1,2,°°+ ©. 
In order that these equations have a solution other than that in which the A; 
are all zero, it is necessary that the determinant of the coefficients be zero. 


On writing the ¢i;(27) simply ¢:;, the determinant is 


1+ o(¢,,—1), oh 2; 
(26) D(o)= 1+ o(¢22—1), = 0. 
This determinant, we see from the fact that the determinant of i; converges 
absolutely, is an integral function of o. Equation (26) is called the funda- 
mental equation associated with the period 27, and does not admit o—0 
as a root. 

As in the characteristic equation in II, the fundamental equation may 
have no finite root, a finite number of finite roots, or a denumerably infinite 
number of finite roots. 

As in the case of constant coefficients, if the fundamental equation has 
no roots, the system (1) has no solution of the form (19), if the fundamental 
equation has a finite number of finite roots, the system (1) has a finite number 
of solutions of the form (19), and finally if the fundamental equation has an 
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infinite number of finite roots, the system (1) has a fundamental set of solu- 
tions, each of the elements of which is of the form (19). 

The discussion of the form of the solutions follows precisely the same lines 
as that in the case of constant coefficients except when o, is an n-fold root of 
(26) and not all the (n—k)th, & > 0, minors of D(c) vanish for o = ap. 

Now let us assume that o = o, is an n-fold root of (26) and not all the 
first minors of D(o) vanish. Then there is only one solution of (1) of the 
form 2; = e%*'y;, where the y; are expressed as in (21). And let us choose 
the notation so that a minor corresponding to the elements of the first column 


is not zero. 
Then we take as a new set of solutions 
(27) Ti, CU" Lij = $ij(t), t= 1,2, °° j=2,3,° ©. 


This set of solutions can be shown to constitute a fundamental set in a pre- 
cisely similar fashion to that in II. Then we make the transformation 


== ettz,, 1,2, ©. 
As above we get 
co 
j=2 


Necessary and sufficient conditions that the 2; be periodic with the period 


are 
24 —2i:(0) = 0, t= 1,2, &. 


On imposing these conditions on (28). we get 


(29) Ay (2) — yi, (0) ] +3 As [pis (2x) — $4; (0)] = 0. 
j=2 

After setting e-*«7 — o/s —1, the equations (29) become 

oo 
(30) = A,(1—o/o,) yi, (0) + jopij + + —1)] = 0, 

kj, 0, +=1,2,--- 

The fundamental equation for the equations (30) is 


* 
(31) (1—o/o,)¥2:(0), 1+ o(¢22.—1), 


Making use of (21) and taking A, —1, equation (31) becomes 


j 
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fo 

sas —(c—1), 

= D(c) = 0. 


So we have 
Y11(0), 


(32) D(c) = (1—o/o,) D, = 1 + 0. 


Since oo, is an n-fold root of D(c) =0, D,(c) has the factor (1— 
o/o,)"*. Since (27) constitutes a fundamental set, any solution can be 
expressed in the form 


foe] 
(33) vi = Byeruty;, + 3 Bids, 4, , 


Now let us make the transformation, corresponding to (19), to get a 
second solution associated with a,. 


(34) ig = (yin + tyi,), 1,2, ©. 


On imposing the condition that x;, shall satisfy the system (1), we find since 
e%*'y;, is a solution, 


co 


From the form of (35) we see that sufficient conditions that the y;, shall be 
periodic with the period 27 are 


(36) — yin (0) =O 3 — 
— ryi,(0) =0. 


On substituting e417 —o,/o, —1, (36) becomes 
co 
(37) —2#(o,— 1) (0) + + Bill + —1] = 0, Aj. 


The condition that these equations be consistent is, since o, 1, 


911(0), °° 
D,(c) = Yo, (0), 1 + 0, (¢2.—1),° =0. 


? 


In (32) we showed that D, vanishes (n—1) times when o=c,. By hypo- 
theses not all the first minors corresponding to the elements of the first column 


i 
38 | 
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are zero for o—o,. Hence we can solve equations (37) for B,, B,,-:-, 
terms of yi,(0). Consequently in this case we get a second solution associ- 


ated with a,, which is of the form (34). 
In a similar manner we can go ahead step by step and get the following 


group of solutions associated with a, 
= 
Gig = [Yin + tyr], 
Gin = [Yin tying 1/(n—1)! 


If o, is a triple root of D(a) = 0, such that all first minors are zero, but not 
all the second minors are zero, the solutions associated with a, are 


iy = Lig = Lig = [Yig + + Yiz) 


All the sub-cases can be treated, as they arise, by the methods given here. 


IV 


Now let us assume that the system of differential equations (1) satisfy 
the hypotheses of the second existence theorem, and in addition the coefficients 
in the power series expansions of the 6;; are separately periodic with the period 
2x. That theorem tells us that the solutions of (1) can be written in the form 


oo 
(39) Vij (t) p*, =1,2,°° ©. 


And we shall take the initial conditions such that 


co 
(40) = Sa (0) = 
whence 
(0) = Cij, (0) = 0, k= 2, 
where the ci; are constants such that c;; = 1, and their determinant is abso- 


lutely convergent. These conditions coupled with the fact that if the deter- 
minant of a set of solutions converges and is not zero when ¢ = 0, it converges 
oo 
and is not zero for every value of ¢ for which 36;; converges, show that the 
4=1 
system (40) constitutes a fundamental set of solutions. 
Now we inquire if we can find solutions of (1) of the form 


(41) erty, 


where the y; are periodic with the period 22, and a is a constant which remains 
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to be determined. After making the transformation (41), the differential 
equations and their solutions become 


oO oo 
(42) 
On imposing the conditions that the y; be periodic with the period 2z, viz., 
yi (2r) —yi(0) = 0, we get 


co co 


After setting — ¢/o —1, the equations (43) become 


co 
(44) O= (2m (2m) wt} — — 1) ] = 0, 


To avoid the trivial case where the A; are all zero, we must set the determinant 
(45) D(o,n) =| [ofa + (x) — c4j(0—1)]| = 0, 


co 
which is a condition on the undetermined constant a. Since the a; ot x0 pk 
k=1 


are the elements of a fundamental set, and ci; = 1, it follows that D(c, ») 
converges absolutely for all finite o’s. If the fundamental equation is satisfied 
when it is also satisfied when where v is any 
integer ; but all distinct solutions of the differential equations can be obtained 
by taking v = 0, since the ratios of the A; are the same for v= 0 as for v= p. 
And for = 0 the equation (45) reduces to 


(46) Do =| (2a) —cij(o —1) | =0. 


When » = 0, we have the case of constant coefficients which we treated 
in II. Now we shall determine the ci; and the a‘ —A- of the solutions as 
we did there. But having the ci; we might also determine the a of the solu- 
tions by means of (46). Since for given initial conditions the solution of the 
differential equations is unique, and since the initial conditions are the same 
in the two cases, the solutions obtained by the former method and the solution 
obtained by means of (46) are the same. 

We recall that the ci; are not all zero if and only if the characteristic 
equation has finite roots. Therefore for every root of the characteristic equa- 
tion there is a corresponding root o of (46); and conversely, for every root 
of (46) there is as corresponding root of the characteristic equation. There- 


H 
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fore if the characteristic equation has no finite roots, then (46) has no finite 
roots, if the former has only a finite number of finite roots, then the latter 
has only a finite number of finite roots, then (46) has an infinite number of 
finite roots which yield independent solutions of the system (1). 

Now since the solution (41) converges for all » sufficiently small, includ- 
ing zero, the failure to find an a for »—0 is not due to clumsy analytic 
methods, but shows that the system (1) has no solution of the form (41). 
From the theory of implicit functions we know that if (46) is satisfied when 
a=a,, then we can solve for a as power series in p, provided » is sufficiently 
small. Therefore the existence of a root of the characteristic equation is both 
necessary and sufficient for the existence of a solution of the system (1) of 
the form (41), and the existence of an infinite number of roots of the charac- 
teristic equation is necessary and sufficient for the existence of a fundamental 
set of solutions of (1), each of the elements of which is of the form (41). 

We shall study in detail the case that the characteristic equation has an 
infinite number of roots; and we shall use only those values of a obtained from 
(45) which for » = 0 reduce to the values of a obtained from the characteristic 
equation. 

When e-?4" — 1, the transformation e°¢™ —o/o—1 can not be made. 
Then as in the case of constant coefficients when a = 0, we can make no gen- 
eral statement about the solution, as the determinant of the A; diverges then. 


Solutions when the a‘ are distinct and a‘ — a‘ =£0 mod V—1. 
The part of (45) which is independent of yu is 


1 — 924,09 
Dy = | — —1) | =| ce; | IH: 
1 — e-2aT 


If (45) were an identity in p, its roots would be the roots of (46), viz., 


a=a‘, Let us assume that we have the general case in which it is not an 


Jj 
identity in pw, and set 


(47) a= + 
Then we get 
(48) D(o, = Dy + (8, 


oo 


e72(a(°) +6,.) 3 e724, (0% 
| ca; | (1 (1 + (a, 
(a0), +5, 


jal 


i 


| 
| 
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where F;, (8, ») is a power series in » and 8, converging for 
| | < co, |» | <p>0. 


Since by hypothesis no two of the a‘ differ by an imaginary integer, the 
expansion of (48) as a power series in 8 and p contains a term in 6 of the 
first degree and no term independent of both » and &. Therefore we know by 
the theory of implicit functions that (48) can be solved uniquely for 4 as a 
power series of the form 


(49) = pPx(n), 


which converges for | » | > 0 but sufficiently small. 

Now we substitute this value of aa‘ + 4, in (44), and get an infinite 
number of linear homogeneous equations for the A; whose determinant con- 
verges and is zero, but the first minors of that determinant are not all zero, 
since by hypothesis, the roots of D, —0 are all distinct and no two differ by 
an imaginary integer. Consequently these equations determine uniquely the 
ratios of the A; as power series in », which converge for p» sufficiently small. 
On substituting these ratios in (42) we have the particular solution yi, 
4=1,2, +++ o, expanded as a power series in wu. Hence we may write it 


oo 
(50) Yyix = (t)p!. 
J=1 
Since the periodicity conditions have been satisfied, 
oo 
yan(t 2x) = (t+ 2x) = 0, 
j= 
for all » sufficiently small and for all real ¢. Therefore 
2m) = 0, 7=0,1,2,°°- 0, 
whence it follows that the y?, j=0,1,2,-° ©, are separately periodic. 
A solution is found in a similar fashion for each a”. 
VI 


Solutions when no two of the a are equal but when a? —a!® 
= 0 mod V—1. 


Suppose that when » = 0 the characteristic equation has two roots such 
that a and a differ by an imaginary integer and that none of the other 
a‘ are congruent to a” mod V—1. Then we see from (45) 
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| Cij | [2 | I] [1 | + (p, ») 
j= 


where as in (47) we have set a—a‘-+ 8,. The term of lowest degree in 6, 
alone is found by expanding the first bracket and turns out to be of the second 
degree. To get the terms in p» alone we suppress those involving 4,, after 
which we get a factor » from each of the first two columns. So we see that 
in general the term of lowest degree in » alone will be in this case of the 
second degree. Hence we have 


1— ear 1 — 
(52) Bam | | E | I] E 
j= 
+ (8;, + &F,(8,, = 0. 


In a similar manner if p of the a‘ are congruent to a mod Y—1, 
then the term of lowest degree in 8, alone is of degree p, and in p alone it is 
of at least the pth degree. 

The problem of the form of the solution of (52) is one of implicit func- 
tions. Writing the first terms explicitly we have 


8? + + + terms of higher degree = 0, 


where Koz) » are constants independent of 8, and On factoring the 
quadratic terms we get 


(53) (8, —d,n) (8, + terms of higher degree 0. 


If d, and d, are distinct, there are two solutions, and these have the 
form * 


(54) = d+ 8:2 = dow + 


where P, and P, are power series which converge for p» sufficiently small. In 
this case the solutions are found as in V. 

But if d, and d, are equal, the character of the solution is in general 
quite different and depends upon the terms of higher degree than the second. 
In general it will be a power series in +VWp. This case we shall consider 
in detail. 

We see from the form of (53) that the expansion of a, as a power series 
in Vp will contain no term in Vp to the first power, but will have the form 


Suppose that this expansion has been obtained from equation (45). 


* Chrystal, “ Algebra,” Vol. 2, pp. 358 ff. 


43 
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Then since a, is not a multiple zero of D, not all the first minors of D are 
zero when a=a,. The ratios of the A; will be determined from (44). If 
yA is a non-vanishing first minor corresponding to an element in the first 
column of D, it follows from the form of (45), remembering that we have 
a(t) = c4jer0#, that solving (44), we get 


pd pa 


where A=A + + 
Aj = ae + 43/2 + AM 

On substituting these series for the A; in (42) we find that the y;, are devel- 

opable as series of the form 


So we see that in general the y;, carry terms in Vz although the term in Vp 
is absent in the expansion for a,. 

However, if all the first minors corresponding to the elements of the first 
column are zero, and if there is a first minor distinct from zero corresponding 
to the elements of the second column, the results are precisely the same. But 
suppose that all the first minors corresponding to the elements of both the 
first and second columns are zero. Then suppose that a first minor corres- 
ponding to an element of the /th column is not zero. Then it follows from 
the form of (45) that when a—ayj, it will carry the factor w?; and let this 
minor be denoted by w?A. Then solving (44) we get 


BA, pra 
A, =— Ax, A; = — Ax, j=3,4,: 
pd 


where A,, A,, 4;, do not in general vanish when »=0. It follows from the 
first two equations that A, must carry p» as a factor, since the A; are finite 
for »=0. Hence in this case the y;, have the same form as before. Simi- 
larly the yi. have the same properties. 

The solutions associated with a‘, a, - + + , are found as in the pre- 
ceding case. If there are several groups of a in which these congruences 
exist the discussion must be made for each one separately. 


Vil 


Solutions when a ts a multiple root. 


Now suppose that two of the a”? are equal and only two, and that there 
are none of the congruences treated in VI. Let us choose the notation so 
that a‘ =a‘. From our work in the theory of linear differential equations 


| 
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in infinitely many variables in constant coefficients and from IV, the solutions 
are of the form 


(54) 3 As Ct pe] or 
(55) = A, 3 yk] + (Cis + +: 3 


After setting o/o—1, the fundamental equation becomes either 
oO 
(56) D(o,n) =| [o +3 2% or 


(51) D(o,#) =| +320? — (o—1)} 
[o{ (Cig + +328 (2a) = 0. 


For »=0 both of these, by the theory of infinite determinants, reduce to 


(58) ey | —o(1— er) [1—o(1— 
g=3 


As before we set aa‘ -+ 6,, and expand as a power series in 6,. Then we 
see that the term of lowest degree in 8, alone is of the second. When the 
determinant is of the form (56) with a‘ —a(, the term of lowest degree 
in » alone is of the second in general. Then we have the same form as in VI. 
But if the determinant is of the form (57), the term in » alone is in general 
of the first degree. In the former case we have a consideration similar to 
that in VI; in the latter case, in general the solutions for 8, are of the form 


(59) = p/?P 
8,2 = P(— p?), 


where P is a power series in Vp and contains a term independent of p. The 
discussion of the special cases is made just as in VI. On substituting these 
expansions for a= a‘- 8, in (44) we solve for the A; as power series in 
Vp. These A; substituted in (42) give yi, and yi, as power series in Vp. 

If p of the a‘® are equal, then for these roots the expansions of D starts 
with 6 as the term of lowest degree in 8, alone, and except in the special cases 
corresponding to those mentioned in the foregoing, the term in yp alone is of 
the first degree. Consequently in general for a = we have 


815 = j= 0, p—l, 


where ¢ is a pth root of unity. 


i 
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VIII 
Solutions when there are equalities and congruences among the a. 


Suppose that two of the a, for example a‘ and a‘) are equal, and that 
a third one, say a‘, differs from a‘” by an imaginary integer. Furthermore 
we shall assume that there are no other equalities or congruences among the 
a. Two cases arise here: (a) the solutions are of the form (54) with 
a‘ =a; (b) the solutions are of the form (55) 

Case (a). In this case we have 


co 
= | Cij | [1 —o(1 — [1 —o(1— 


In setting aa + 8,, we find that the term in D of lowest degree in 6, 
alone is of the third. To get the term in D of lowest degree in p alone we se* 
8, = 0. Then it becomes evident at once that each of the first columns carry 
p as a factor, while the remaining ones do not. Consequently the term of 
lowest degree in » is of the third degree at least. Furthermore since the first 
three columns vanish when » —8,—0 there are no terms of lower degree 
than the third in 6, and w. Hence in general we see D of the form 


The problem is now one of implicit functions. The details of the special cases 
must be treated as they arise. However, we make the general statement that 
since the roots of the cubic terms of (60) set equal to zero are in general 
distinct, it follows from the theory of implicit functions that the three values 
of 6, are in general expansible in integral powers of p. 

Case (b). In this case we have 


co 
Dy = [1—o(1 — [1 —o(1 — — 0, 
j=4 


On introducing 8, as before, we find that the term of lowest degree in 8, alone 
is of the third degree. But when the terms involving p are retained in D, only 
the first and third columns vanish when » 8,0, and consequently the 
expansion of D will contain a term in p? alone. Furthermore, since the first 
and third columns vanish for » = 8,—0, there will be no terms of degree 
lower than the second in » and 8,. Hence in general D has the form 


(61) D=8 + + + 0, 


In the general case in which y,, and y,, are not zero, there is one solution in 


integral powers of » and two in powers of Vz. 
When the roots a‘ have higher multiplicities and more congruences 


among them, we make a similar discussion. 
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Periodic Orbits on a Surface of Revolution. * 


By DANIEL BUCHANAN. 


§1. Introduction. The object of this paper is to determine the periodic 
orbits described by a particle which moves, subject to gravity, on a smooth 
surface of revolution, the axis of which is vertical. Let us denote the equation 
of the surface by 


2? + y?—2pa-+ (22+ y%, 2) =0, (1) 


where p is a positive constant, e an arbitrary parameter, and ¢ a power series 
in + y? and z, converging for 2? -+ and |2| sufficiently small. There 
will be no loss of generality if we suppose that the constant term of ¢ is zero, 
otherwise it could be eliminated by a linear substitution for z, that is, by a 
translation of the zy-plane along the z-axis. Since the equation (1) vanishes 
with x? + y? and z, it may be solved, by the theory of implicit functions, for 
a? + y? as a power series in 2 converging for | z| sufficiently small. Hence 
the equation (1) may be expressed as 


F(a, y,2) =a? + y? 2ef(z) =0, (2) 


where f is a power series in 2 converging for | z| sufficiently small. 

For e= 0 the equation (2) represents a vertical paraboloid of revolution. 
The generating parabola has its axis coinciding with the positive z-axis and 
its semi-latus-rectum equal to p. Periodic orbits are first constructed when 
the particle moves on this paraboloid. Then the analytic continuation of 
these orbits is made with respect to « and in this way the orbits for the more 
general surface (2) are determined. It was for this reason that the parameter 
e was introduced in (2). 

The form of the equation of the surface (2) was suggested by a somewhat 
similar equation used by Poincaré in his memoir, Sur les Lignes Géodésiques 
des Surfaces Convexes: Trans. Am. Math. Soc., vol. vi (July, 1905). Geodesic 
lines play the same réle in Poincaré’s memoir as periodic orbits do in this 
paper. 

It will be readily seen that when «0 the problem here considered is 
somewhat analogous to the well-known problem of the spherical pendulum.t 


* Presented to the American Mathematical Society, Sept. 5, 1918. 

+ For a complete discussion of the spherical pendulum, including the horizuntal 
as well as the vertical motion, see Moulton’s Periodic Orbits, chap. III, also Rendiconti 
cel Circolo Matematico di Palermo, Vol. XXXII (1911), pp. 338-365. 
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§2. The Differential Equations. If the particle is of unit mass and 
moves without friction, then the differential equations of motion are 


y"=—Y, 


where X, Y, and Z are the normal reactions due to the surface, and g is the 
acceleration due to gravity. Since the surface is assumed to be smooth, the 
normal reactions at any point are proportional to the direction cosines of the 
normal at that point, and therefore the differential equations become 


= Y = = (3) 
2’ = Z—g = p+ fz) — 9, 
where A is a factor of proportionality. 
These equations admit the vis viva integral 
a? + 2? 2g(c—2z), (4) 


where c is the constant of integration. 

Since the equation of constraint (2) does not contain ¢ explicitly, the 
factor A can be obtained by differentiating F'( x, y, z) twice with respect to ¢ 
and eliminating 2’, y”, and 2” from the result by means of the differential 


equations (3). Then on making use of the relations 
+ y? — Ref, + = — 2” + 29g(c—2z), 
from (2) and (4), respectively, we find 


__ —9(p— + 2 (1 — efzz) (5) 
p? + 2pz—2(f+ 


The part of 2A which is independent of « is 


2A 


+ 22) 


Xo 


(A). Prrtopic OrBITs ON A PARABOLOID OF REVOLUTION. 


§3. Proof of Existence of Periodic Orbits. For «= 0 the surface (2) 
becomes 
x? + y? — 2pz = 0, (7) 


and the differential equations (3) become 


=) =Ay, 2” =— J, (8) 


H 
H 

i 
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where A, is defined in (6). The last equation of (8) is independent of the 
first two and admits the integral 


— + Cy 
(9), 


where ¢, is the constant of integration. 

If the vertical motion is to be periodic, z cannot increase indefinitely and 
therefore 2’ must vanish for some value of t, {= ¢), say. Suppose z= 2%, at 
t=—t,. Then since 2 —0 at t= 1%), it follows that the constant of integra- 
tion in (9) has the value 


Cy == — 49% (¢ — %). 


Hence the integral (9) becomes 


(% — 2) (4 —c-+2). (10) 
It is readily seen from (10) that 2’ vanishes for z= 4%, and also for 
Z=C— %. 

As p is assumed to be positive, no part of the paraboloid (7) will lie 
below the zy-plane. Then z and 2 cannot be negative and therefore, for real 
initial conditions, c cannot be less than z,. Three cases arise according to the 
values assigned to c. They are: 


Case I. < 
Case IT. 
Case III. c> 


Case I. Suppose 4% Sc < 2% or c—%<%. Then # is real so long 
as 2 lies in the interval 2, >z2>c—2. If the particle is started in the plane 
2=2Z, it cannot remain in this plane for zz, +40 does not satisfy the 
z-equation of (8). Moreover, it cannot move above this plane since z < %, 
hence it must move below this plane. Consequently z decreases or 2’ is nega- 
tive, and the negative sign must be taken in (10). The particle continues to 
fall until it reaches the plane z= c— 2 where the velocity again vanishes. 
As z=c—42, is not a solution of the z-equation in (8) and as z cannot be 
less than c — 2, it must increase, or the positive sign must be taken in (10). 
The particle then rises until the plane z=, is reached where the velocity 
again vanishes and the radical in (10) changes sign. Hence the particle 
oscillates between the two horizontal planes z = 2 and z= c — %, the latter 
plane being the lower plane. 


4 
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Case II. If c—22 or c—% = 4%, then 2 would be real only when 
2==%, and would therefore be identically zero. In this case the particle 
would revolve in the horizontal plane z=, the centrifugal force and the 
reaction of the surface being just sufficient to overcome gravity. 


Case III. If c>2z or c—%> % it can be shown by an argument 
similar to that used in discussing Case I that the particle oscillates between 
the same two planes z = 2 and 2 = c — 4, as in Case I, but the latter plane 
is the higher plane. 

Thus if ¢ is positive and not less than 2, the vertical motion of the 
particle is periodic, and therefore the last equation of (8) must admit a 
periodic solution. 

It could also be shown by the usual analytic existence proof that the 
last equation of (8) admits a periodic solution. It is not necessary, however, 
to establish the existence of a periodic solution by either method for, by 
Macmillan’s theorem, quoted below, it is shown that if the formal construc- 
tion of a periodic solution can be made, then this solution will converge under 
suitable conditions. Macmillan’s theorem is as follows: * 

If 


d 
Tn} pt) (t=1, (11) 


is a system of differential equations in which the fi are expansible as power 
sertes in , Un, and p, vanishing for +: = p= 0, with coeffi- 
cients which are uniform, continuous, and periodic functions of t with the 
period 2jx: and if the fi converge for 0 St S 2jx, when | < pj, |u| <r 
then the solutions 2;(t) are expansible as power series in p, or any fractional 
power of 4, which converge for all t in the interval 0S t S 2jx provided | p| 
is sufficiently small. If the constants of integration can be determined at 
each step so as to make the solution formally periodic with the period 2jx, then 
the solutions so determined will be periodic and converge for all finite values 
of t provided | | is sufficiently small. 


§4. The Vertical Motion in terms of Elliptic Integrals. The integral 
(10) may be written 


We desire to express this integral in terms of elliptic integrals. The neces- 


* Trans. Am. Math. Soc., Vol. XIII, No. 2, pp. 146-158. 
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sary substitution is found to be 


Vmn(1—xw) 
where 
m= p+ 2%, n= 2z, 
Vm—vVn 


Since p, c, and % are all positive, then m is also positive. If z< p/2+c 
then n is likewise positive, and therefore x is real and less than 1 numerically. 
When the above transformation is made in (10’) we obtain 


2V mn (1—xw) dw 
=2V9(t—t), 
Vm aid Gow V9( ) 

which simplifies to 


2V mn 2k 1 — w’? 
Vm + Vn 2 3 1 K2 2V g(t to), 


where 
dw 
V (1 — w?) (1 — 


an elliptic integral of the first kind, and 


dw 
(1—x2w?) V (1— w?) (1 — 


an elliptic integral of the third kind. Thus ¢ can be found for any given 
value of w or z. For most purposes, however, it is more convenient to have z 
expressed as a function of ¢. This will be discussed in the next section. 

The half-period, 7/2, of a complete oscillation can be obtained by 
integrating (10’) between the limits z, and c—z,. The corresponding limits 
for w are —1 and +1. Hence 


— mn dw 
2 T /2) = 


+ dw 
J. (1 — x2w?) V (1 — w?) (1 — 
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Since PES = 0, then 
4\/mn dw 


d 


§ 5. Construction of a Periodic Solution for the Vertical Motion. Let 


2=v-+¢/2, (12) 
where v is the new dependent variable. Then the last equation of (8) becomes 

+2’ 


where k=p-+c. This equation can be expanded in the form 


= — (1/k) (4gv + [1 — 20/k + 40? /k? 
(13) 


which will converge for |v| <|&/2|, or z< p/2+c. Hence the particle 
must not rise above the plane z = p/2 +c, being the plane at the distance ¢ 
above the horizontal plane passing through the focus of the generating para- 
bola of (7). 

Now let 


t—t=4Vk/g (14+ 8) 7, (14) 

where y and 8 are arbitrary parameters. When (14) is substituted in (13), 
the factor y can be cancelled off in the resulting equation, and we obtain 

ib = —[(1 +8) w+ piv?) J, (15) 


where » = y/k, and the dots denote derivation with respect to r. 

The terms of (15) which are independent of » and 6 are i + w—0, 
and the period of the solution is 2x in r. Now an existence proof, based on 
Poincaré’s extension to Cauchy’s theorem, would show that a periodic solution 


of (15) exists having the form, 


(16) 


ms 


fo 
w= wip, 
j=0 
where the w; are periodic functions of + with the period 2, and the 4; are 
constants; and that this solution converges for all values of + in the interval 
0 S17 S 2z, provided | » | is sufficiently small. Instead of making this rather 
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simple existence proof, we shall assume the form of solution (16), make the 
formal construction, and then show that Macmillan’s theorem will apply to 
prove the convergence of the solution. 
Since 
+ c/2, (17) 


and since 2 = 0 at {= then at r—0 (or we may choose, 
w(0) =a, w(0) = 0, 


where a is a real constant. When these initial conditions are imposed on 
(16), we have 


W (0) =a, = 0, j= 1," = 0, j= 0,1,-° 00. (18) 


Let (16) be substituted in (15) and let the resulting equation be 
denoted by (15’). ‘This equation is an identity in » and we may equate the 
coefficients of the same powers of p, thus obtaining differential equations which 
define the various w;. These equations are to be integrated, and the constants 
of integration and the various 8; are to be chosen so that the solutions will be 
periodic and satisfy the initial conditions (18). 

Step 0: Coefficient of p®. The terms of (15’) that are independent of 
p give the equation i, + w,—0, and its solution which satisfies (18) is 
Wy = COS T. 

Step 1: Coefficient of ». The differential equation obtained from the 


terms in p is 

i, + Wy= — 8, — = — a8, cos + (a?/2)(1 + 3 cos 27).(19) 
In order that w, shall be periodic the coefficient of cos r in the right member 
of (19) must be zero, otherwise the solution would contain the non-periodic 
or Poisson term rsinr. On putting 6,0, the general solution of (19) 


becomes 
w, = A, cost + B, sin r + (a?/2) (1— cos 2r). 


From the initial conditions (18) it follows that the constants of integration 
A, and B, are both zero, and the desired solution at this step becomes 


8, = 0, w, = (a*/2) (1 — cos 27). 
Step 2: Coefficient of w?. The differential equation at this step is 


We + We 4wyW, 4w,° Wo" 2WyW, 
= — a(8, + %a?) cost — cos 38r, 
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and the desired solution is 
8, = — ya’, w, = cos + cos 3r). 


The remaining steps of the integration are entirely similar to the pre- 
ceding step and by an induction to the general term it can readily be shown 
that the process can be carried on indefinitely. It will be observed that 4; is 
zero if j is odd, and that each w; carries the factor a/+* and is a sum of cosines 
of multiples of + having the opposite parity of 7. The highest multiple of r 
in w; is (j+1). 

So far as the computation has been made, the solution of (15) is 


243 
w=acost + /2) (1 — cos 2r)— (cos cos (20) 
§ = — +- 


The convergence of this solution will now be considered. In order to 
put the differential equation (15) in the same form as (11) of Macmillan’s 
theorem, we let 

w= W,, dW,/dr¥ = W., 
then (15) becomes 


adW,/dr = f,(W,, W.3u,7) = W,, 
aW,/dr = f.(W,, Wo; m7) =—[(1 +5) W, + 


The f, and f, are obviously expansible as power series in W,, W., and p», and 
vanish for W, = W, =p=—0. The coefficients are constants and therefore 
satisfy the condition that they shall be uniform, continuous, and periodic 
functions of r. Further, f, and f, converge for W,, W., and p sufficiently 
small numerically. Hence equations (15) satisfy the conditions of Macmil- 
lan’s theorem and consequently the solution (20) converges for | »| suffi- 
ciently small. 

When (20) is substituted for w in (17) we obtain 
16 
As » and @ are both arbitrary and as they occur only in products as indicated, 
we may suppress either without loss of generality. Let us suppose a—1. 
Then the periodic solution for the vertical motion of the particle becomes 


k[p cos + (p?/2) (1 — cos 27) — (73/16) (cos r — cos 3r) 


— + ( 


272 
0/2 + k[nacos + (1 — cos 2r) (cos ]. 


| 
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In the expansion for z, the multiples of + in the coefficient of »/ have the same 
parity as 7 and therefore 

= ¢/2—kp=c/2—y. 
Further 

2(0) 
Hence the particle oscillates between the planes z = c/2 + y and z= c¢/2—y. 
The parameter y is therefore a scale factor denoting the amplitude of the 
oscillation on either side of the plane z= c/2. 


§6. The Horizontal Motion. When the solution (21) has been substi- 
tuted in the first two equations of (8) and the transformation 


has been made, we obtain 


and the same equation in y, where 
6, = — cost, 06, = — 14(1 — cos 2r), 


and the remaining 6; are likewise sums of cosines of multiples of + having 


the same parity as 7, the highest multiple being 7. 

Equation (22) is similar to the one first discussed by Hill* in 1877 
in his celebrated memoir on the motion of the lunar perigee. A very complete 
list of references to the literature of the differential equations of this type is 
given by Baker on page 134 of his memoir “ On Certain Linear Differential 


Equations of Astronomical Interest.” + 
Three cases arise in the solution of (22) depending upon the values of 


(k/4p).~ They are 
Case I. 4Vk/p 0 and Vk/p not an integer. 
Case II. 4Vk/p0 and Vk/p an integer. 
Case III. 4Vk/p=0. 


Since k = p-ce and since p and ¢ are both positive in the physical 


problem under consideration, it follows that 4Vk/pF 0 and Case III need 
not be considered. 


*The Collected Works of G. W. Hill, Vol. I, pp. 243-270; Acta Mathematica, 
Vol. VIII, pp. 1-36. 

t Phil. Trans. of the Royal Society of London, Series A, Vol. 216, pp. 129-186. 

t'Compare also Moulton’s Periodic Orbits, Chap. ITI, § 52. 
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Case I.—This may be regarded as the general case. 


The form of the solution of (22), first given by Floquet,* is 


etary, (23) 


where a and wu are power series in p, the former having constant coefficients, 
the latter periodic coefficients with the period 27 in r. When (23) is substi- 
tuted in (22) we obtain 


ti + 2a — a?u + k/plt + +--+ Ju=0. (24) 
Now let 


(25) 
be substituted in (24). We obtain a differential equation (24’), which is an 
identity in w. Since the right number is zero, then the coefficient of each 
power of » must also be zero. On solving the various differential equations 
thus obtained, we determine the uw; and aj, the uw; as periodic functions of +r 
having the period 27, and the aj; as constants so determined that the w; shall 
be periodic. 
Since the solution (23) is later multiplied by an arbitrary constant, see 
equations (39), we may choose u(0) 1, from which it follows that 


=1, uj(0) =0, @. (26) 
Step 0: Coefficient of w°. From the terms of (24’) that are inde- 
pendent of », we have the differential equation 
tin +iVk/p 
and its solution is 


Up = + (27) 


where a, and by are the constants of integration. Since Vk/p is not an 


integer in this case, the term e~*V#/pr does not have the period 2x and we put 
b= 0. From (26) it follows that a1 and the desired solution at this 
step becomes uy = 1. 


Step 1: Coefficient of wu. On equating to zero the terms in p in (24’) 
we obtain 


ti, iV /p tt, = Vi /p a +k/p cosr. _ (28) 


* Annales de Vv Ecole Normale Supérieure, 1883-4. 
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The complementary function is 


U, =a, -+ 


a, and b, being arbitrary constants. 

It is well known in the theory of differential equations that the presence 
of any term in the right member which has evactly the same period as any 
term of the complementary function will yield Poisson terms in the particular 
integral, that is, terms containing + outside of trigonometric or exponential 
symbols. Since the constant part of the right member of (28) has exactly 
the same period as the constant in the complementary function, then non- 
periodic terms will arise in the particular integral unless a, —0. On putting 
a, =0 the complete solution of (28) becomes 


From the periodicity and initial conditions it follows that 


k 


b, = 0, 4, 
p—k 


and the desired solution at this step becomes 


uy = [1 — + sin tr], a, = 0. (30) 


Step 2: Coefficient of w?. The differential equation at this step is 


°p(p—k) 
k? 
+ (cos r— $cos 2r + sin 27) — (k/p) cos2r , 


lig + k/p =a,Vk/p T (31) 


and the solution which satisfies the periodicity and initial conditions is 


4p(p—k)?(k — 4p) (p —k)? (cost —iVk/p sin r) 
— 4p) cos 27 + 1Vk/p 4(p—k)(k —4p) sin 27, 
—— (5k—p 


So far as the computation has been carried out it is observed that the 


5Y 
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aj are real constants and that the u; consist of sums of cosines and i times 
sines of multiples of 7, the highest multiple being 7. It will now be shown 
by induction that these properties hold in general. 

Let us suppose that a,,° °° , ans, Uo, Ui, °° * » Un have been com- 
puted and that the various a; are real and that the uj, 7—=1,°--, n—1, 
have the form 

j 
Uj= cos Ir + sin Ir), (33) 

1=0 
where the _ and B‘” are real constants. We wish to show that an 
can be determined as a real constant and that uw, has the same form as (33) 


when j =n. 
Step n: Coefficient of w™. The differential equation at this step is 


The only undetermined constant which enters (34) is a» and it is written 
explicitly in so far as it occurs. The function U, is linear in u,,°-* , Un_y- 
The terms of Un which arise from 2iaii in (24) have the form 


and are cosines and ¢ times sines of multiples of r, the highest multiple being 
nm—1. The terms which arise from a?u in (24) have the form 


Jar am l,m = 0, n—1, 


where ( ) denotes 1 if 1m and 2 if 14m. These terms have the same 
form as the preceding terms. When / =m = 0, the term —a,?up will cancel 
off with the term (4/4p)u» arising from the last factor of (24) and therefore 
the highest multiple of 7 in the sines and cosines is m»—1. The last factor of 
(24) gives, in addition to the term (k/4p)un just considered, terms of 
the type 


(k/p) 


These terms are also of the same form as (33) but the highest multiple of r 
which they yield is n. Therefore 


n 
Un= [a™ cos lr + sin Ir], (35) 
1=0 


where a™ and b™ are real constants. 
Let us now determine the periodic solution of (34). In order that the 


i 

i 
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solution shall be periodic, the right member of (34) must contain no constant 


terms. Hence 
dn = Vp/ka™, 


a real constant. The complete solution of (34) then becomes 
Un = An + + Un, 


where @, and b, are the constants of integration, and the particular integral 
Un has the same form as (35), except that it has no constant term. From 
the periodicity and the initial conditions we have 


bn = 0, a= — Un(0), 


and the desired solution for un is the same as when j=. Hence the prop- 
erties of aj and wu; already stated hold in general. 

A second solution of (22) could be constructed in an entirely similar 
way but it is not necessary to make the construction as this solution can be 


obtained directly from the former solution. 
Let the solution of (24) already obtained be denoted by u(r, +7). 


Then one solution of (22) is 


ettar y (7, (36) 


Since the differential equation (22) is independent of 1, a change in the sign 
of iin (36) will still give a solution, viz., 


2 = e-tar u(r, —i). (3) 


Thus a second solution can be obtained by changing the sign of 7 in the first 


solution. 
The determinant formed by the two solutions (36) and (37) together 


with their derivatives with respect to 7 is a constant,* and its value can be 
computed most readily when r—0. This determinant is 


D=— 2[ia + &(0, + =—iVk/p + terms in yp, (38) 


which is different from zero for » 0 and therefore remains different from 
zero for | »| sufficiently small. The two solutions (36) and (37) therefore 
constitute a fundamental set and the most general solution of (22), as well 
as of the similar equation in y, is 


* Moulton’s Periodic Orbits, § 18. 


| 
i 
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x= A,etaty™ + 
y= + B, 


== u(r, +4) = 1+ —cost + iVk/psint]u 


4 —22pk? + 
4p(p—k)?(k — 4p) (p—k)? 
2pk? — 4p — 
4p(p—k) (k — 4p) 
u® (+ 1) =u (—i4), 


where A,, A,, B,, and B, are arbitrary constants. 

By using Macmillan’s theorem, quoted in § 3, it can readily be shown 
that the above solutions converge for all r in the interval 0 = r S 2z provided 
| «| is sufficiently small. 


(cost —iVk/p sin 7) 


COs 27 


Case II. Suppose Vk/p =v, an integer. Sincek=p-+c and p and 
c are both positive, then y >1. Two sub-cases arise depending upon the 
value of vy. They are: 


Sub-case I. v=2, 
Sub-case II. v2. 


Sub-case I. The construction proceeds as in Case I until equation (27) 
is reached, and this becomes 
Uy = Ay + dye-?*7, (40) 


where both terms have the period 27. From the initial condition u,.(0) =1 
it follows that b, = 1— ap, and the solution for wu, becomes 


Uy = My + (1 (41) 
with the constant a) remaining arbitrary at this step. 
Step 1: Coefficient of ». When we equate to zero the coefficient of p 


in (24), that is in equation (24) after (25) has been substituted and Vk/p 
replaced by 2, we obtain 


ti, + 2a,a) — 2a, (1 — ay) + 2a,e*7 
e-tr + 2(1 —a,) #7, (42) 


! 
| 
| 
| | 
J 


f 
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The complementary function of this equation is 
Uy, = a, + 


Since terms of exactly the same period as those of the complementary func- 
tion occur in the right member of (42), viz., constants and terms in e-?7, 
the particular integral will contain non-periodic terms unless the constants 
and the coefficient of e*7 are put equal to zero. Since we seek a periodic 
solution we therefore put 


0, 2a,(1—a) = 0. (43) 


These equations are satisfied by a, 0, a arbitrary, but a, must be different 
from 1 or this case would be the same as Case I. The complete solution of 
(42) then becomes 


U, = a, + — eit + — % (1 —a,) (44) 


Since u,(0) = 0, then 
b,=— (4% +4). 


At this step a and a, still remain undetermined. 


Step 2: Coefficient of w?. The terms of (24) which contain the factor 
py? give the equation 


lig + = + + 2] + [— 2a, (1 — ay) — + 1% 
+ 2a,e't — — 
— 2(a, + %)e*r — 1% (1 — a, 


On equating to zero the constants and the coefficient of e*‘7 in the right 
member, as at the preceding step, we have 


Ray — + 2—0, (45) 
— 2a,(1— a) — %a + 1% =0. 

These equations are satisfied by 
ag + % = YaV2l?, (46) 


where the upper signs are to be taken together, also the lower. When equa- 
tions (46) are satisfied, the general solution for uw, will be periodic, being 


Uz = + — + na e747 — 


+ (a, + %)e** + 
Since u,(0) = 0 then 


| 
| 
| 
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b, = — a, + — *%e- 
The undetermined constants at this step are a, and 4. 
Step 3: Coefficient of »®. The differential equation at this step is 


— U9, + u,(4 — 4008 27) + 4u, cos 


If we consider only the part of U; which involves the arbitrary constants that 
are determined at this step we have 


Us = 2agM + a, (2a, + %) +d, 
+ [—2(1 —ap)a; + @, (2a, — %) + de (47) 


+ terms in e*!7, e747, e-tir, 


where d, and d, are known constants. To make wu, periodic we equate to zero 
the constants and the coefficient of e*‘7 in (47) giving 


Raz + (2a, + %) +d,—0, (48) 
— 2(1—ay)as + a, (2a, — %) + d, = 0. 
The determinant of the coefficients of a, and a, in (48) is 


where the upper sign here is to be taken with the upper signs in (46), simi- 
larly for the lower signs. Since this determinant D, is different from zero, 
equations (48) can be solved for a, and a,, and the solutions are unique for 
either set of values in (46). When (48) is satisfied, the solution for u, will 
therefore be periodic. Two constants of integration, a, and b,, say, will 
arise from the complementary function. One constant, b,, say, will be deter- 
mined from the initial condition u,(0) —0, while a, and a, remain undeter- 
mined at this step. 

The succeeding steps of the integration are entirely similar to step 3 
just considered. So far as the constants of integration are concerned one of 
them, b; say, is determined by the initial condition uj(0) —0 at the step 
where it arises. The other constant a; of the step j, is not determined until 
the step 7 +2 is reached, where, on equating to zero the constants and the 
coefficient of e-*7 in the right member of the differential equation in w;,,, 
two linear equations in a;,, and a; are obtained and the determinant of their 
coefficients is the same as D, in (49). Hence a;,, and aj can be uniquely 
determined provided one set of values is taken in (46). 


ti, + — Ug = — Piagtiy — -+ Ragtty 
| 
D, = + %V21%, 
| 


PO 
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It is thus evident that two solutions of (24’’) can be obtained in the 
foregoing construction, according as the upper or lower signs are taken in 
(46). The corresponding solutions of (22) are 


+(%+ 
and 


If we proceed as in Case I to obtain two solutions by changing the sign 
of 7 in the preceding solutions it would appear that two additional solutions 
could be obtained, but this is impossible since the differential equation (22) 
is only of the second order and admits of but two solutions. This apparent 
difficulty is overcome if the factor e‘7 is multiplied into the part of the solu- 
tions contained in the square brackets [ |]. Thus 


= et — 1,217) et + 
+ (% + + € Jaber (50) 

+ 


It is obvious, therefore, that the two solutions v, and 2, differ only in the sign 
of or of 

If the exponentials in the square brackets of (50) are expressed in 
trigonometric form, these solutions become 


= Lo == tary (2), 

a= YoV217 p?+- 3 (51) 
u™ —=costr—t¥V 217 sinr- + ( | 

(+ i) =u (—i). 


The terms represented by wu“ and w® differ not only in ¢ but also in the sign 


of ¥217. They are power series in » with sums of cosines and iV/217 times 
sines of multiples of +, the highest multiple in the coefficient of p/ being 
j+1. 

The determinant of the two solutions in (51) together with their deriv- 
atives is a constant, as in Case I, and its value at r—0 is 


D, = —2[ia + (0)] =i % ( 


| 


| 
4 
4 
2 
| a 
| 
i 
| 
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which is different from zero for » = 0, and therefore remains different from 
zero for |» | sufficiently small. Hence the two solutions (51) constitute a 
fundamental set, and the most general solutions of (22) and the similar 
equation in y are 

A,etary™ + Anetary®, y= Bietary™ + (52) 


where a, wu, and wu are defined in (51), and A,, A,, B,, B, are the constants 
of integration. 

Sub-case II. When v is an integer different from 2, the construction 
is the same as in the preceding sub-case until equations (43) are reached. 
The equations analogous to (43) are 


= 0, va,(1— ad) = 0. (53) 


These are satisfied by a, = 0, a, arbitrary, but a, must be different from 1 as 


in the preceding sub-case. 
The equations which correspond to (45) are 


(54) 
v(5v7—1)\ _ 
(1—a,) { \ == (), J 
Since a, 1 in Case II, then the solutions of (54) are 
__ v(5v? — 1) 
The equations analogous to (48) are 
v? (5y? — 1) 
a, { + \ == (), — vas + v6 = 0, (56) 


where 6° is the constant part of 0, in equation (22). The solutions of these 
equations are 
a, = 0, a, = v8. 

The constants a, and a, are determined in the same way at the next 
step as a, and a, were found in the preceding step. They have the same 
coefficient in the equations analogous to (56) as a, and a, have in (56), and 
similarly for the succeeding steps. It will be found that 


aj = 0, j= 0,1,- 


and that the remaining a; are, in general, different from zero. 


& 
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So far as the computation has been made we find 


Uo =e 


v? 


-ip-lr 


v+1 v—2 v+2 


When these terms are substituted in (25) and the result in (23) we obtain 
one solution of (22). Let it be denoted by 


== etary), (58) 
Another solution can be obtained by changing the sign of 7 in (58), thus 


etary, (2), (59) 


where u® 4) =u (—1). 


The determinant of these two solutions and their first derivatives is 


which is different from zero for | »| sufficiently small. The solutions (58) 
and (59) therefore constitute a fundamental set, and the general solutions 
of (22) and the corresponding equation in y are of the same form as (39) 


or (52), viz., 
z= A,etary™ +. y Betaty™ + Bie-tary), (60) 


where a and uw have the values in (58). If the exponentials in (58) are 
expressed in trigonometric form then 


uP — [AH + BY cos(v + + sin(v + 1) p/, (61) 


j=0 1=0 


where A‘, BYP, and C‘ are real constants. In particular 


AP = 0, j= 0,1, 2, BY CY, j=0,1, 2, l= 0,1, 2, 
v? 


2  2(v—1)’ 3 2(v+1)’ 
5 


| 

n | 
r 
| 

| 
| 
| 
| 
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BO =x + 
0 “(F—1)? + 2(v—1)(v?—1) ’ 
2(v-+1)(%—1)’ 32(v— 2)’ 


B32 + 2)" 


From (57) it follows that 
AP —=0, 7=0,1,--° , 


The remaining A‘ are different from zero, in general. 

If the factors e*/®)r and e-»/2)r are taken with u®™ and w, respect- 
ively, in (60), then u® and uw will have cosines and sines of (v/2 +1)r 
instead of (v + 1)r as in (61). 


§ 7%. The Arbitrary Constants. Each of the solutions for x, y, and z 
contains two arbitrary constants, and the values of these constants for the 
physical problem will now be considered. 

In this discussion no discrimination will be made between the two cases 
of the preceding section. It was for this reason that the same notation was 
chosen for the solutions (39), (52), and (60). 

Besides satisfying the differential equations (8), the solutions for z, y, 
and z must also satisfy (7), viz., 


a? + = 0, (7) 


and their derivatives with respect to ¢ must satisfy the vis viva integral (4). 
When the transformation 


t— tp = 4Vk/g(1+8)r, 8=— 
is made in (4), this integral becomes 
e+ y+ = (k/2)(1+8) (c—z). (62) 


Since any horizontal section of the surface (7) is a circle with centre 
on the z-axis, the zy-axes may be rotated about the z-axis without changing the 
form of (7), and we may therefore suppose that the particle lies in the xz-plane 
at the initial time, or y= 0 at r—=0. Since (0) = (0) —1, then it 
follows that B, + B,=—0. 

As 2(0) =0, the particle must be initially projected in a horizontal plane. 
Since any horizontal section of the surface of constraint is a circle with centre 
on the z-axis and since u(0) = 0, then it follows that When this 
condition is imposed on the solution for z, we have 


| 
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[ia + [A,—A,] = 0. 
The factor [ta + wi (0) ] is different from zero for | » | sufficiently small and 


therefore 
A, = A,. 


Hence the solutions for x and y become 
A[etary) +4 e-tary?)], y Bletary® — etary], (63) 


where the now unnecessary subscripts on the constants A and B have been 


dropped. 
On putting r= 0 in (63) and (21) we get 
= 2A, =0, 
y(0) =90, y(0) = | (64) 
2(0) =c/2 + kp, 4(0) = 0. 


When these values are substituted in (7) and (62) we find 


A= +4Vp(c+ = +4Vp(c-+ 2y), 
__ , VE +8) 2), 
4fta + (0) J 4fa— wi™ (0) | 


Thus A is real and B is purely imaginary. Hence the solutions of equations 
(8), in terms of r, are 


tVk(1 +8) (C—2y) retary) — etary (2) 


+H 


2=c/2+k[p cos + — cos 27) (65) 
=— 
n=y/k, k=p+e. 


The terms a, wu, and uw are defined in (39), (52), and (60) according to 


the value of Vk/p. The double signs in (65) depend upon the octants of 
space into which the particle is initially projected. 

If the exponentials in (65) are expressed in trigonometric form, the 
solutions of (22) for x and y become 


t= Vp(ce+ 2y) [{ },cosar—{ },s8in ar], 


; 
| 
| 
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where { }, and { }, are power series in » which have different values for the 
different cases of §6. For Case I, 


co j 


{ =1 + Las Cos Ir |p, 


j=1 I= 


co j 
{ & sin 
j=1 I=1 


where aj; and bj; are the coefficients of cos /r and sin Ir, respectively, in w;, 
equations (30), (32), and (33). 


For Case II, 
j 
{ 3 [AM 4 BY cos(v + 
1=0 
co j 
{ }o=3 [OY cin(v + 1) 
j=1 


The constants A‘?, BY’, and CY’ have the same values as in (61) for sub- 
case II, and for sub-case I they are, in so far as the computation has been 


carried out, 
AM = 0, BO 1, CO 217. 


There still remain in the solutions (65) two constants which have not 
been determined. They are c and y (=p). Since c is the constant of 
integration arising in the vis viva integral, its value depends upon the initial 
velocity of the particle. As 7’ —2’ —0 at = the initial velocity may be 
denoted by yf. Then from the vts viva integral (4) it follows that 


=g(c— 2y), or C— Ry = 
Since, further, c + 2y = 2z, then 
C= + y=4(%— 9/7/29) - 


Thus the constants c and y are functions of the initial velocity and the initial 
height. 
Let us now return to the three cases of §3. Since 2 —c/2-+ y these 
cases become: 
Case. OC ySc/2, 
Case II. y=0, 
Case III. y<0, |y|Sc/2. 


In Case I and Case III the particle oscillates between the horizontal planes 
2—=c/2-+y and z=c/2—y. The two orbits are geometrically the same 
but the motion in the one orbit is half a period ahead of the motion in the 
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other. If y= c/2, then the initial velocity is zero and as there is no lateral 
projection the particle will move in the vertical parabola 2?—2pz. If 
y=—c/2, then 2 = 0, and the particle is projected from the lowest point 
with the initial velocity y’7—= V2gc. It will therefore move in the vertical 
parabola y? = 2pz which is dynamically the same orbit as 77 2pz. The 
orbits when y = +c/2 correspond therefore to the simple pendulum. 

In Case II when y= 0, then z—c/2, a constant, and the differential 


equations of motion become 
a” + (g/p)a=0, y” + (g/p)y = 0, 2” =0. 
Their periodic solutions which satisfy the initial conditions are 
V pe cosV g/p t, you V pe sin V g/p t, 2z=c/2. 


In this case the particle moves in a circle the plane of which is parallel to the 
xy-plane and at a distance c/2 above it. 

In the first paragraph of § 5 it is stated that the inequality z << p/2-+¢ 
must hold in order that the expression in (13) will converge. Now the 
maximum value of z is 


a=c/2+]|y|, 
and therefore the above inequality becomes 


or ly|<3(p+¢). 
Since | y | must not exceed c/2 in order that the initial velocity shall be real, 


it follows that the inequality z < p/2-+c will always be satisfied for real 
initial conditions. 


(B). Prrropic OrBits oN A SURFACE OF REVOLUTION. 


§ 8. The Method of Solution—Let us next consider the construction 
of the periodic orbits described on the more general surface of revolution 
represented by the equation 


F(a, y, 2) = 2? + y? — 2pz + = 0. (2) 


The differential equations of motion have already been found, equations 
(3) and (5). The method of constructing the periodic solution of these 
equations is to first make the analytic continuation with respect to « of the 
periodic solution for the vertical motion obtained in (A) where e=0, and 
then substitute this solution for z in the first two equations in (3). We thus 


| 
| 
| | 
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obtain two differential equations having periodic coefficients somewhat similar 
to (22). 


§9. The Equation of Varvation—Let us substitute in the last equation 
of (3) 
Z2=2-+ (66) 
where Z denotes the solution obtained in (21), and é is a function of + which 
vanishes with «. We obtain 


(67) 


where the undefined terms have the following properties: 

(1). The functions ©, and ®, are periodic functions of 7, but we shall 
show that it is not necessary to know the explicit values of these functions in 
order to solve (67). 

(2). The functions Z,, Z,,-- - , Z; denote power series in é having 
coefficients which are power series in » with sums of cosines of multiples of + 


in their coefficients. These functions also contain additional terms in é and &, 
the former being multiplied by power series in » with sums of sines of mul- 
tiples of + in the coefficients, the latter by similar series except that they 


contain cosines. The function Z, contains no linear terms in é or é. 
If we neglect the right member of (67) we obtain 


é+0,é+ 0,¢=0, (68) 


which is called the equation of variation. The generating solution is z =z, 
or the expression for z in (21). 

Now it has been shown by Poincaré * that if the generating solution 
contains an arbitrary constant which does not occur in the original differential 
equations of motion, viz., equations (3), then a solution of the equations of 
variation can be obtained by differentiating the generating solution with 
respect to this constant. Three constants occur in (21), viz., c, w, and f,, the 
latter entering implicitly through +r. The constant ¢ occurs in 2A, equation 
(5), and therefore enters the differential equations (3). The two remaining 
constants, » and ¢, do not occur in (3) and therefore may be used in applying 
Poincaré’s theorem. Each constant yields a solution of (68) and since the 
differential equation is only of the second order, both its solutions may be 
obtained from the generating solution. Hence it is not necessary to know 
the coefficients ©, and @, in (68) in order to solve the differential equation. 


* Les Méthodes Nouvelles de la Mécanique Celeste, Vol. I, chap. tv. Loc. cit. 


— 
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Consider, first, the constant ¢. Then one solution of (68) is 


where 
S(r) =sin r— psin 27 — Agp? (sins —3sin3r) +--+. (70) 


Since this solution is later multiplied by an arbitrary constant, see equation 
(76), we may drop the constant factor of S(r) in (69) and take 


8(r) (71) 
as the solution. 
Considering the constant p, we obtain as the second solution of (68) 


02 Or 08 


(72) 


where the parentheses ( ) denote differentiation in so far as » occurs explicitly 
inZ Now 


=kC(r) =k[cos + — cos 27) 
2 
(cost —cos 3r) +++: (73) 
—} 1 -+ power series in p?]. 


Therefore the solution (72) becomes 


&=k[C(r) + 7KS(r)], 
where 
K = p?[— 1 + power series in p?]. 


As in the previous solution we may drop the constant factor & and take 
+ 7KS8(r), (74) 


as the second solution of (68). 

The two solutions (71) and (74) constitute a fundamental set, since the 
determinant of these two solutions together with their derivatives is different 
from zero for | »| sufficiently small, being 


A =—1-+ power series in p. (75) 
Hence the general solution of (68) is 
é=N,S(r) + N.[C(r) + rK8(r)], (76) 


where N, and N, are arbitrary constants. 


ice 
| 
| | 
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From the way in which S(r) and C(r) were derived, it is readily seen 
that the coefficients of sin(j + 1)z7 and cos(j-+1)r in the coefficients of p/ 
in S(r) and C(r), respectively, are equal. 


§ 10. Construction of the Solution for the Vertical Motion. 
Let 
(77) 


where the €; are to be periodic with the period 2z in r. 
Since the initial time was chosen for «0 so that z(0) =0, it may 
now be chosen so that 2(0)—0 when «+0. Now z—2+6, and as 


=0, it follows that (0) =0. Hence 
&(0)=0, (78) 


Let (77) be substituted in (67) and let the resulting equation be 
denoted by (67’). It is an identity in e, and on equating the coefficients of 
the same powers of e we obtain sets of differential equations which can be 
integrated and the various constants of integration can be chosen, as we shall 
show, so as to satisfy the periodicity and initial conditions. 

Two types of series occur in this integration and they are similar to 
S(r) and C(r) in (71) and (%3), respectively. These series will be 
denoted by 


S;(r), SP (r), 8i(r), (r), j=1, 
or 
C;(r7), CP (r), C;(r), CH (7), j=1, 


according as they are similar to S(r) or C(r), respectively. 
Coefficients of «. When the coefficients of ¢ to the first power are equated 


in (67’) we obtain 
é, + @,é, + ©,é, =O (7). (79) 


The complementary function of this equation is the same as the solution of 
the equation of variation, viz., 


& + +7 KS(r)], (80) 


where n™ and n™ are the arbitrary constants. By employing the method 
of the variation of parameters to find the complete solution, we have 


(r) +H [O(r) +7K S(r)] =0, 


(81) 


— 


BucHaNaNn: Periodic Orbits on a Surface of Revolution. %3 


The determinant of the coefficients of %? and 7% in the above equations is 
the same as (75), and since it is different from zero equations (79) can be 
solved for and Thus 
WO = — (1/A)C™ (7) [C(r) +7 K S(r)], (82) 
(1/A)C™ (7) S(r). 
Since the coefficients of the same power of » in C™(r) and C(r) are sums 


of cosines of the same multiples of 7, the product C®(r)C(r) will yield, in 
addition to periodic terms, a constant 


where p are real constants. Then the integration of (82) gives 


n® —N® — [pr py +8 (r) + 
No (r), 


2 


(83) 


where VN and N® are the constants of integration. Since the coefficients 
of the sines and cosines of the highest multiples of + in the coefficients of p/ 
in S(r) and C(r), respectively, are equal, then it follows that S®(r) and 
C™(r) have the same property. When (83) is substituted in (80) the 
complete solution of (79) then becomes 


&, = NYS (r) + NO [O(r) + — + (1/H) (7), 


where C,(r) contains no terms independent of yw. In order that é, shall be 
periodic V® must be given the value 


vp, 1 
(power series in p), (84) 
and from the initial conditions (78) it follows that 
N® = 0. 
Hence 
i= (1/p) C, (r), | (85) 


where 0,(r), like C,(r), contains no terms independent of p. 
Coefficients of «®. Equating the coefficients of <? in (6%) gives the 
differential equation 


+ = (1/u)O™ (86) 
Except for the factor 1/p in the right member, this equation is similar to 
(79). The general solution is obtained in the same way as at the preceding 
step and is found to be 
6 


hl 
= 
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= NS (r) +N + +772 + (2), 


where N® and N‘ are the constants of integration, and p, is a power 
series in w. In order to satisfy the periodicity and initial conditions we 
must put 


P2 1 
(power series in »), =0. 


Hence 
= (1/n?) C2 


where C’,(7), contains no terms independent of p. 
The remaining steps of the integration can be carried on in the same 
way and by an induction to the general term it can be shown that 


n= (1/n")Cn(r), 


where C,(r) as at the previous steps, contains no terms independent of p. 
On substituting the solutions for the various é; in (77) we obtain 


(81) 


By applying Macmillan’s theorem, quoted in § 3, it is found that the solution 
(87) converges for | «| sufficiently small. 

When (87) is substituted in (66), the solution for the vertical motion 
becomes 


Cy (7) (88) 
Since ¢ and p are both arbitrary we may put « = py and therefore (88) becomes 
2=Zz+3 C; (7) pi. (89) 
This solution converges for | »| and | p| sufficiently small. 


§10. The Horizontal Motion. When the first two equations of (3) 
are transformed by the substitution 
t—t—=4v 


already used in (A), and the value of z obtained in (89) is substituted in 2a, 
these differential equations become 


(90) 


t+ [bo + dip + + ]r=0, \ 
9 + + bop? +° ly=9, 
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where each ¢; is a power series in » with sums of cosines of multiples of + 
in the coefficients. The function ¢, has the same value as the coefficient of 
g in (22). 

Now let pop. Then the coefficients in (90) can be rearranged as 
power series in » and the differential equations take the same form as (22). 
Hence the solutions are 


L, + L, e*Prv,, y = M, e*Prv, + M, (91) 


where L,, L,, M,, and M, are the constants of integration, and f, v,, and v, 
are similar in form to a, u,, and w,, respectively, of (39), (52), or (60) 


according to the value of Vk/p. 
When the initial values +(0) —y(0) —0 are imposed on (91), it 
follows that 


[iB + 0,(0)] [Z,—L,] =9, M,+M,=—0. 
Since 18 + 1,(0)£0 for | »| sufficiently small, then 
L, = L, = L, say, 
and, from the second equation, 
M, =— M,=M, say. 
Hence the solutions (91) become 
x= L[etBry, + e-Bry,], y= — e-*Bry,]. 


The constants Z and M can be determined as in § 7, and it is found that L is 
real while M is purely imaginary. Suppose M=—iN. Then 


co j j 
2L[cos Br{1 +3 & piarcos — sin Br{S sin I7}], 
j=1 1=0 j=1 
y = 2N[sin Br{1 & pia cos Ir} — sin Br{S sin Ir} ], 
j=1 1=0 j=1 
where a‘? and 6‘ are real constants. These equations together with equation 
(89) are the solutions of the equations of motion of the problem under 


consideration. 
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